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Abstract. We study Neumann type boundary value problems for nonlocal 
equations related to Levy processes. Since these equations are nonlocal, Neu- 
mann type problems can be obtained in many ways, depending on the kind of 
"reflection" we impose on the outside jumps. To focus on the new phenom- 
enas and ideas, we consider different models of reflection and rather general 
non-symmetric Levy measures, but only simple linear equations in half-space 
domains. We derive the Neumann/reflection problems through a truncation 
procedure on the Levy measure, and then we develop a viscosity solution the- 
ory which includes comparison, existence, and some regularity results. For 
problems involving fractional Laplacian type operators like e.g. (— A) Q / 2 , 
we prove that solutions of all our nonlocal Neumann problems converge as 
a — > 2~ to the solution of a classical Neumann problem. The reflection mod- 
els we consider include cases where the underlying Levy processes are reflected, 
projected, and/or censored upon exiting the domain. 



1. Introduction 

In the classical probabilistic approach to elliptic and parabolic partial differ- 
ential equations via Feynman-Kac formulas, it is well-known that Neumann type 
boundary conditions are associated to stochastic processes having a reflection on 
the boundary. We refer the reader to the book of Freidlin [9j for an introduction 
and to Lions and Sznitman 17J for general results. A key result in this direction 
is roughly speaking the following: for a PDE with Neumann or oblique boundary 
conditions, there is a unique underlying reflection process, and any consistent ap- 
proximation will converge to it in the limit (see [17] and Barles & Lions [5]). At 
least in the case of normal reflections, this result is strongly connected to the study 
of the Skorohod problem and relies on the underlying stochastic processes being 
continuous. 

The starting point of this article was to address the same question, but now 
for jump diffusion processes related to partial integro-differential equations (PIDEs 
in short). What is a reflection for such processes, and is a PIDE with Neumann 
boundary conditions naturally connected to a reflection process? It turns out that 
the situation is more complicated in this setting, at least the questions have to be 
reformulated in a slightly different way. In this article we address these questions 
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through an analytical PIDE approach where we keep in mind the idea of having a 
reflection process but without defining it precisely or even proving its existence. 

For jump processes which are discontinuous and may exit a domain without 
touching its boundary, it turns out that there are many ways to define a "reflec- 
tion" or a "process with a reflection" . This remains true even if we restrict ourselves 
to a mechanism which is connected to a Neumann boundary condition (see below). 
But because of the way the PIDE and the process are related, defining a reflection 
on the boundary will change the equation inside the domain. This is a new nonlocal 
phenomenon which is not encountered in the case of continuous processes and PDEs. 

PIDE with Neumann-type boundary condition - In order to simplify the 
presentation of paper and focus on the main new ideas and phenomenas, we consider 
different models of reflections and rather general non-symmetric Levy measures, but 
only for problems involving linear equations set in simple domains. The cases we 
will consider already have interesting features and difficulties. To be precise, we 
consider half space domains ft :— {(xi, ...,xn ) = (x',x N ) e R N : x N > 0} and 
simple linear Neumann type problems that we write as 



(1.1) 



i(x)-I[u](x)-f(x) = in O, 

dx 



9u = in dn, 



or sometimes as 



J F(x,u,I[u}) = in Q, 

|-#^ = o in an, 

K ox N ' 

where F(x, r,l) = r — I — f(x) and 

I[u](x) = lim / [u(x + rj(x, z)) — u(x)] dfi(z). 



b— J-0+ 



b<\z\ 



We will assume that / <E Cf,(fi), i.e. / is bounded and continuous, that /i is a 
nonnegative Radon measure satisfying 



(1.2) / \z\ 2 A 1 dn{z) < oo, 
and that 

(1.3) x + r](x, z) £ Q for all x G Q , rj(x, z) = z if x + z G 17. 



Note that I[u] is a principal value (P.V.) integral, and that (1.2) is the most 
general integrability assumption satisfied by the Levy measure associated to any 
Levy process pQ. When r)(x, z) = z, then I[u] is the generator of a stochastic 
process which can jump from x € f2 to x + z with a certain probability, see e.g. 



[U [HI US] . Assumption (1.3) is a type of reflection condition preventing the jump- 
process from leaving the domain: nothing happens and T)(x, z) = z if x + z G f2, 
while if x + z (jt then a "reflection" is performed in order to move the particle 
back to a point P(x, z) = x + r](x, z) inside fl. Note that we have to check at some 
point that the reflection is consistent with a Neumann boundary condition. 
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The main examples of r\ arc the following model cases, where we use the notation 

x = (x' , xn) <= R^ -1 x M + , rj(x, z) = (rj(x, z)' , rj(x, z)n), etc.: 



(a) i](x,z) 

(b) r)(x, z) 

(c) rj(x,z)-- 

(d) 77(1,2;) 



2 if xn + zn > 
if not 

[ z if xn + zjy > 
I zf^i if not 

z if xat + zat > 

(z', — xn) if not 



if Xn + z N > 



[censored] 

[fleas on the window] 
[normal projection] 
[mirror reflection] 



(z! . — 2xn — zn) if not 
for all x € fi and z ^ 0. The different reflections are depicted in the figure below. 

Censored rj a Mirror reflection r\d 

Flea on the window ?^ Normal projection rj c 



x + r)i(x, z) 
f 



X + 7]i(,X, Z) 



V 

x + z 



Dashed lines are jnmps 



We will discuss later whether the naively proposed "reflections" are realized by a 
concrete Markov process, i.e. if they correspond to the generator of such a process. 

Main results - From an analytical (PIDE) point of view, we first have to give a 



sense to problem (1.1) and relate it to an homogeneous Neumann boundary value 



problem. This is done in Sections [2] and [3] The first part is classical: to take into 
account singular Levy measures, we split the integral operator in two and write the 
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equation in a more convenient way. Here classical arguments in viscosity solution 
theory are used, see e.g. |3] and references therein. Viscosity solution theory is also 
used to give a good definition of the Neumann boundary conditions. 

If fj, is a nice bounded measure, then the problem ( 1.1 ) can be solved easily with- 
out caring much about the Neumann boundary condition. Moreover, the solutions 
will be uniformly bounded by ||/||oo- Intuitively (1.1) carries the information that 
the particles remain in Q since they can only jump inside Q. This mass conservation 
is an other way to understand that we are dealing with a (homogeneous) Neumann 
type of boundary condition. 

When fi then is a singular measure, we can approximate it by a sequence of 
bounded measures (/z n )„, consider the associated (uniformly bounded!) solutions 
(u n ) n , and wonder what the limiting problem is if fi n converges to /z in a suitable 
sense. This is the way we choose to make sense of both the definition of problem 
(1.1 ) and the associated notion of (viscosity) solutions. We point out here that the 
"real" Neumann boundary condition arises only if the measure is singular enough. 
In the other cases, either the process will never reach the boundary as in the 
censored case for a-stable process with a < 1 (see e.g. 0), or the equation will 
hold up to the boundary. 

The natural next step is then to prove uniqueness results for all the above models 
and equations. In this paper different types of proofs are given depending on the 
singularity of the measure and the structure of the "reflection" mechanism at the 
boundary. These results are given in Sections g] - 16| The first case we treat is 
when the jump function r\ enjoys a contraction property in the normal direction. 
This covers all the non-censored cases listed above - models (b)-(d) . Had we had a 
contraction in all directions, then the usual viscosity solution doubling of variables 
argument would work. Here we have to modify that argument to take into account 
the special role of the normal direction. 

The second case we consider is the censored case (a) when the singularity of the 
measure is not too strong. By this we mean typically a stable process with Levy 
measure with density -J^ ~ | z |"+a for a €E (0,1). We construct an approximate 
subsolution which blows up at the boundary and this allows us to derive a compar- 
ison result by a penalization procedure. Such a construction is related to the fact 
that the process does not reach the boundary in this case, see e.g. [B]. 

The last case is the censored case (a) when the singularity is strong, i.e. when 
a G [1,2). This case requires much more work because no blow-up subsolutions ex- 
ist here. In fact, when a £ [1, 2), the censored process does reach the boundary (see 
e.g. [5]). We first prove that the Neumann boundary condition is already encoded 
in the equation under the additional assumption that the solution is /3-H61der con- 
tinuous at the boundary for some /3 > a — 1. Then we prove a comparison result for 
sub/super solutions with this Holder regularity at the boundary. The proof is then 
similar to the proof in the a < 1 case, except that the special subsolutions in this 
case are bounded and only penalize the boundary when the sub/super solutions are 
Holder continuous there. Finally, we construct solutions in this class. In dimen- 
sion N = 1, we use and prove that any bounded uniformly continuous solution is 
Holder continuous provided fi satisfies some additional integrability condition. In 
higher dimensions, we use and prove a similar result under additional regularity 
assumptions on / in the tangential directions. 
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Finally, in Section [7] we show that all the proposed nonlocal models converge to 
the same local Neumann problem when the Levy measure approches the "local" case 
a = 2. More precisely, we consider Levy measures fi a with densities (2 — a) upif-jL ; 
where g is a nonnegative bounded function which is C 1 in a neighborhood of and 
satisfies g(0) > 0. In this case we prove that whatever nonlocal Neumann model 
we use, the solutions u a converge as a — > 2 to the unique viscosity solution of the 
same limit problem, namely 

-aAu — b ■ Du + u = f in 17 , 

(L4) < du n ■ flo 

— = mail, 

on 

where a := g(0) 1 -^-^ and b := Dg(0) 1 -^-^ . 

Related work - One of the first papers on the subject was Menaldi and Robin 
[15] . In that paper stochastic differential equations for reflection problems are 
solved in the case of diffusion processes with inside jumps, i.e. when there are no 
jumps outside the domain. They use the method of "penalization on the domain" 
inspired by Lions, Menaldi and Sznitman [18]. 

In model (a) (the censored case), any outside jump of the underlying process 
is cancelled (censored) and the process is restarted (resurrected) at the origin of 
that jump. We refer to e.g. [TO] EH US] for more details on such processes. The 
process can be constructed using the Ikeda-Nagasawa-Watanabe piecing together 
procedure [6] QUI [16] , as Hunt processes associated to some Dirichlet forms [6] [13] , 
or via the Feynman-Kac transform involving the killing measure [6 . In particular, 
the underlying processes in this paper arc related to the censored stable processes 
of Bogdan et al. [5] and the reflected a-stable process of Guan and Ma [TB . But 
note that we essentially only construct the generators and not yet the processes 
themselves. On the technical side, we use viscosity solution methods, while [H [13] 
use the theory of Dirichlet forms and potential theory. Our assumptions are also 
different, e.g. with our arguments we treat more general measures and we have the 
potential to treat non-linear problems, while jS] 113) e.g. treat much more general 
domains. 

Let us also mention that the "natural" Neumann boundary condition for the 
reflected ct-stable process of Guan and Ma [TJ] is slighly different from the one 
we consider here. They claim that the boundary condition arising through the 
variational formulation and Green type of formulas is 

limt 2 - a ^{x + te N ) =0. 
t->o dx N 

This formula allows the normal derivative to explode less rapidly than |a;jv| a_2 - 
In model (b) outside jumps are stopped where the jump path hits the boundary, 
and then the process is restarted there. Model (c) is close to the approach of Lions- 
Sznitman in |17j . and here outside jumps are immediately projected to the boundary 
along the normal direction. This type of models will be thoroughly investigated in 
the forthcoming paper [3] by three of the authors, but this time in the setting of fully 
non-linear equations set in general domains. Note that model (b) and (c) coincide 
in one dimension, i.e. when N = 1. Finally, in model (d), outside jumps are mirror 
reflected at the boundary. This is intuitively the natural way of understanding a 
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"reflection", but the model may be problematic to handle in general domains due 
to the possibility of multiple reflections. E.g. it is not clear to us if it correspond 
to an underlying Markov process in a general domain. 

To the best of our knowledge, processes with generators of the form (b)-(d) have 
not been considered before. E.g. the works of Stroock jT9] and Taira [20l |21] seem 
not to cover our cases because their integrodifferential operators involve measures 
and jump vectors rj that are more regular than ours. Moreover in these works, it is 
the measure and not rj that prevents the process to jump outside f2. 

In the case of symmetric a-stable processes (a subordinated Brownian motion), 
our formulation follows after a "reflection" on the boundary. So such processes 
can be constructed from a Brownian motion by first subordinating the process 
and then reflecting it. Another possible way to construct a "reflected" process in 
this case would be to reflect the Brownian motion first and then subordinate the 
reflected process. A related approach is described e.g. in Hsu [T3] where pure jump 
processes like Levy processes are connected via the Dirichlet-Neumann operator 
to the trace at the boundary of a Reflected Brownian Motion in one extra space 
dimension £1 x K + . An analytic PIDE version of this approach is introduced by 
Caffarelli and Silvestre in [7] in order to obtain Harnack inequalities for solutions of 
integrodifferential equations, and then these ideas have been used by many authors 
since. 

Finally we mention the more classical work of Garroni and Menaldi , where a 
large class of uniformly elliptic integro differential equations are considered. There 
are two main differences with our work: (i) In the principle part of all equations 
is a local non-degenerate 2nd order term. This allows the non-local terms to be 
controlled by local terms (the solution and its 1st and 2nd derivatives) via interpo- 
lation inequalities, and the local W 2 ' p and C 2 ' a theories can therefore be extended 
to this nonlocal case. In our paper, it is the non-local terms that are the principal 
terms, and interpolation is not available. In addition, most of our results can be 
extended to degenerate problems, (ii) In [111 Dirichlet type problems are consid- 
ered, and the authors have to assume extra decay properties of the jump vector rj 
near the boundary, conditions that are not satisfied in our Neumann models. 



2. Assumptions and Definition of solutions 



In this section we state the assumptions on the problem (1.1 1 that we will use 



in the rest of the paper, give the definition of solutions for ( 1.1 1, and show that the 



quantities in this dcfintion is well-defined under our assumptions. 

In this paper we let Du(x) and D 2 u{x) denote the gradient and hessian matrix 
of a function u at x. We also define P(x, z) = x + rj(x, z), and then we can state 
our assumptions as follows: 

(%) feC b (Q). 

(H M ) The measure fj, is the sum of two nonnegative Radon measures and 



where c is either or 1, /1, is a symmetric measure satisfying (1.2) and 



\z\ d[i* 



and 



\z\<l 



(1 A \z\) dfi# < 00. 
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(EL) Neuman problem: P(x, z)n — xn + T)(x, z)n > for any x, z and 

rj{x, z) = z for any xn + Zn > 0. 
(H*) At most linear growth of the jumps: there exists c, ; > such that 

\i](x, z)\ < c n \z\ for any x, z. 
(H?) Antisymmetry with respect to the z'-variables: for any i = 1, . . . , (N — 1), 

r)(x, <Jiz)i = -r)(x, z)i where <j % z = . . . , -z i} z N ) . 
(H^) Weak continuity condition: 

r)(y, z) rj(x, z) /i-a.e. as y x. 
(H^) Continuity in the x'-variable: 

\tj(x', s, z) — r](y , s, z) | < C|z| |x — y | for any x , y' , z and any s > 0. 
(Hjj) Non-censored cases: Contraction in the iV-th direction: 

\P(x, z) N - P(y, z) N \ < \x N - y N \. 
(H^) Censored case: For all z 7^ and a; € ft, 



r)(x,z) 



z ii xn + zn > 0, 
otherwise. 



Remark 2.1. Assumption (H M ) means that we can decompose [i into a sum of a sin- 
gular symmetric Levy measure and a not so singular Levy measure. Symmetric here 
means that J RN xi z ) d^i = for any odd /x-integrable function \- This assumption 
covers the stable, the tempered stable, and the larger class of self-decomposable 
processes in R. , cf. chapter 1.2 in [I]. In all these cases the Levy measures satisfy 

d\i g(z) 



dz \z\ N + a 



for a G(0,2), 



and (H M ) holds with c = if a £ (0, 1), while if a € [1, 2) and g is Lipschitz in 
-Bi(O) and bounded, then (H^) holds with c = 1. In the last case we may take 

= j^l and ^# = 9 M^m for Mz) := min g{yl 

dz \ Z \ N + a dz \ Z \ N + a ^ ' |2,| = |*| 

and note that h is symmetric and g — h is nonncgative. More generally, we can 
consider measures where 4^ = g(z)-^ t and /Lt* is symmetric. 

Remark 2.2. The cases (a), (b), (c), and (d) mentioned in the introduction satisfy 
assumptions (Hi) for i = 0, 1, 2, 3, 4, where in fact (ELJ holds with C = 0. Assump- 
tion (H„) holds except in case (a), and case (a) is equivalent to (H®). Note that 77 
is continuous in x for z 7^ in (b), (c), and (d), while in (a), r\ is continuous except 
on the codimension 1 hypersurface {zn — xn}- 

Now we will define generalized solutions in the viscosity sense, and to do that 
we need the following notation: 

I[tf>] =I 5 [ct>} + I s [<fi}, 

where 

I s [<p} = (t>{x + rj(x, z)) - <)>{x) dn(z) . 

J\z\>6 
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The I s -term is well-defined for any bounded function <j>. For the /,5-term there are 
two cases, depending on whether c = or 1 in (H M ). If c = 0, a Taylor expansion 
shows that Ig[(/)](x) is well-defined for <f> G C 1 and If c= 1, and the measure 

/i is very singular, we add and subtract a compensator and write 



for 



I 5 [cj>}(x)=I s [<t>}(x)+P.V. / D<f>(x)-t](x,z)dn(z), 

J\z\<8 



Is[<f>}(x) := / <f>(x + 77(0:, z)) — <f>(x) — D<j>(x) ■ r](x, z) dfi(z). 

J\z\<S 



By the C -regularity of 0, these two terms will be well-defined - see Lemma 2.1 



below. Note that this results is non-trivial because the compensator is not well 
defined in general! 

Definition 2.1. Assume that (HJ, (H*,) for i = 0, 1, 2 hold. 

(i) A bounded use function u is a viscosity subsolution to if, for any test- 

function <fi G C 2 (IR Ar ) and maximum point x of u — <j> in -B(a;, c^<5) n O, 

F(z,u(x),/ 5 [0] + / <5 [u]) < ifa;GO and 

either F(x,u(x),I s [4>] + I s [u}) < if x e 50 and c = 0, 

-(a;) < if x G 90 and c = 1. 



(ii) A bounded lsc function u is a viscosity supersolution to (1.1) if, for any test- 
function <fi G C 2 (IR Ar ) and minimum point x of v — cf> in -B(x, c^c)) ("1 O, 

F(af,u(af),/i[</>] +/ 5 H) > if x G and 

either F(x, ^[^ + J 5 ^]) > if x G dCl and c = 0, 

(x)>0 if x G 90 and c = 1. 



9a; at 

(hi) A viscosity solution is both a sub- and a supersolution. 

Remark 2.3. The constant c, ; is defined in (Hi). If u and </> are smooth and a; is a 
maximum point of u — (f> over B(x, c v S) D O, then by (Hi), 

u(x) — <j>(x) > u(a; + r/(x, z)) — <p{x + 7]{x, z)) for all \z\ < S. 

If we rewrite this inequality and integrate, we find formally that 7,s[t(](x) < I$[<f>](x). 
Lemma |2.1 1 below makes this computation rigorous. From this inequality it is easy 



to prove that classical (sub) solutions of ( 1.1 ) are viscosity (sub)solutions. Moreover, 
smooth viscosity (sub)solutions are classical (sub)solutions (simply take <p — u). 

Remark 2.4. In general to pose boundary value problems in the viscosity sense, 
one requires that either the minimum of the equation and the boundary condition 
is nonpositive or the maximum of the equation and the boundary condition is 
nonnegative. Here this is not the case and for a natural reason. If the measure is 
very singular (and c = 1) then the equation cannot hold on the boundary and the 
inequality holds just for the boundary condition. In the c = case, on the contrary, 
the equation will hold up to the boundary and the boundary condition can not be 
imposed in general. In other words, we only end up with a Neumann boundary 
condition if c = 1, i.e. the measure has a "strong" singular part /x*. In this case 
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the intensity of small jumps is so strong that the jump-reflection mechanisms, e.g. 
as in (a) - (cQ, are not enough to prevent the process from "diffusing" onto the 
boundary, and we need to add a reflection process at the boundary to keep the 
process inside (just as in the case of Brownian motion). We also note that the 
symmetry of /z* is a natural condition in order to obtain Neumann and not oblique 
derivative boundary conditions, cf. Lemma |3.3| and proof. 



We now prove that Is[4>] is well-defined for <f> £ C 2 . 

Lemma 2.1. Assume (H M ) and (HI) for i = 0, 1,2, and let x € Q, <f> £ C 2 , and 
5 > 0. Then I$[(f>](x) is well-defined since 



I s [<j>](x) = I s [<l>]{x)+P.V. / D<f>(x)'ti(x,z)dii(z), 

J\z\<8 

and the compensator term satisfies 



P.V. / D(j){x) ■r 1 {x,z)dp(z) 

J\z\<5 

Dcf>(x) ■ rj(x, z) d[i#(z) + c / Dcf>(x) ■ rj(x,z) dfj,*(z). 

\z\<6 Jx N <\z\<8 

Moreover, there is R = R(x, 77) > such that 

Is[<f>](x) = S (l)\\<p\\ C 2(B R (x)) as S^O. 

In the following, we often drop the P.V. notation for such integrals since they 
may be expressed in terms of converging integrals. Note that the integral over 
{xjv < \z\ < 6} need not vanish since this regions exceeds the boundary and hence 
rj(x, z) will not be odd there. 

To prove Lemma |2.1[ we need the following result. 



Lemma 2.2. Assume (H M ) and (HI) fori = 0,1,2, and let xn > 0, p £ (0,xn), 
and v £ R N be a fixed vector. 

(i) For r £ (0, p), j r< \ z \ <p v ■ rj{x, z) dp(z) = J r<N<p v ■ r)(x, z) dp # {z), and 



P.V. / v ■ T](x, z) dp(z) — I v ■ n(x, z) dp#(z) . 
J\z\<p J\z\<r> 

{ii) Forr £ (0,1), J r<z \ <s v' ■ rj(x, z)' dfx(z) = f r< \ z \ <s v' ■ rj(x, z)' dp, # (z), and 



P.V. / v' ■ i](x, z)' dp(z) — / v r ■ rj(x, z)' dp,#(z) , 

J\z\<8 ' J\z\<8 

(Hi) For r £ (0, 1), 

rj(x, z)n dp*(z) > I (zn — xn) dp*(z) > . 



r<\z\<8 



r<| 2 |<5 
N >3: N 



Proof, (i) If \z\ < p < xjv, then r)(x, z) = z by (H ). Hence r\ is odd with respect 
to the z variable, and the integral with respect to the symmetric part c/z* is zero. 
Passing to the limit as r —> and using the integrability of p# along with (H^) 
finishes the proof of (i). 
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(ii) Let <j be the rotation o~(z(zn) = (— z(zn)- Since ji„ is symmetric, 



v' ■ rj(x,az)' d/j,*(z) = / v' ■ rj(x, z)' d[i*(z). 

<|z|<5 Jr<\z\<8 

Other hand, since rj(x, —z', zjv)' = —T](x, z', zjv)' by (H 2 ), the above integral is zero. 
The result on the principal value is obtained as in the first case, after letting r — > 0. 
(iii) Let us decompose 



r<|z|<<5 



r](x, z)n d[i*{z) = I r)(x,z)N dfi*(z) 



r<{z{<5 



T](x,z)n dfu(z) + / rj(x, z) N dn*(z) . 

r-<|2|<4 J r<\z\<6 

Z N> X N Z N<~ X N 

The integral over — xn < zn < xn vanishes since T](y, z) — z in this region and /z* 
is symmetric. By (H?) we have that T)(x,z)n > —xn if zn < —xn and rj(x,z) = 
zn > xn if zjv > xn- Hence by symmetry of /i*, 



rj(x, z)n d(i*(z) > / (zn — xn) dfi*(z) > 0, 



<\*{<t 

iJV>a=iV 



'r<|z|<<5 

and the proof is complete. □ 



Proof of Lemma \2.1\ The expression for Is is obtained by adding and substracting 
the compensator term. The first integral in this expression is well-defined since 
the integrand is smooth and bounded by the function \\z\ 2 max B (j. m |-D 2 0|, for 
R = max y6S ( 0;( 5) \r)(x, y)\, which is an /i-integrable function over B(0, S). Moreover, 
/o<|z|<<5 l z l 2 — as 5 — > since \z\ 2 is /x-integrable near 0. 

In the compensator term, the integral with respect to n# exists by (H*), while 
the integral with respect to /x* over B(0,xn) vanishes by Lemma [2.2| -(z). Since \z\ 
is integrable near the origin for /z^, this term is \D<p(x)\os(l) as S — > 0. □ 

3. Derivation of the boundary value problem - PIDE approach 

In this section we derive the boundary value problems from approximate prob- 
lems involving a sequence of bounded positive Radon measures /x fc = l{|z|>i/fc}M 
converging to [i. Assume (H M ) and let ^ = l{| z |>i/ fe }M# and ^ = l{| z |>i/fc}M*) 
it then easily follows that 

(Hi) Jra^y" \z\ A 1 d^(z) = j\z\ A 1 d^(z), 
(Hi) Jirn^ j \z\ 2 A 1 d^{z) = [ \z\ 2 A 1 d^(z), 



(Hi) lira / \z\ A 1 = oo. 

fc— >+oo J 

The approximation problem we consider is then given by 
(3.1) u(x) - I m [u}(x) = f(x) in H, 

where, for <j> 6 Cb(f2), 

W$0) = / 0(^ + ??(a;,z)) - 4>{x) dn k (z). 
J\z\>a 
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Since the measures [i k are bounded, this equation holds in a classical, pointwise 
sense. Moreover, it is well-posed in C{,(f2) and the solutions U}. are uniformly 
bounded in k: 

Lemma 3.1. Assume (H/), (H M ), (H°), and (H^). 



(a,) i*br every k, there is a unique pointwise solution Uk of (3.1) in C(,(f2). 

(7),) //ufc and ffc are pointwise sub- and supersolutions of (3.1), then Uk < Vk in tt. 

(c) If Uk is a pointwise solution of (3.1), then \\uk\\L x (n) < ||/ Hi/WO)- 

Proof, (a) Let T : Ct,(il) — > Cb(f2) be the operator defined by 

Tu := u - e(u - 1^ [u] - f) , 

where e<(l + 2||/i fe ||i) _1 and H^lli is the total (finite!) mass of the measure fi k . 
Then T is a contraction in the Banach space Cf,(f2) since 

\\Tu-Tv\loc < (l-e)||w-w|| 00 + 2e|| A i fe || 1 || u -w|| co 



< 



(l-e(l + 2\\^\\ 1 ))\\u-v\ 



KCWWu-vWoo, 

and C(k) < 1. Hence there exists a unique Uk G Cf,(il) such that Tuk — Uk, which 
is equivalent to ( |3.1| . 

(b) If supq(m — v) is attained at a point x € fi, then by the equation and the 
easy fact that I Mfc [4>] < at a maximum point of 0, 

sup(u — v) = u(x) — v(x) < [u — v](x) < 0. 
n 

The general case follows after a standard penalization argument. 

(c) Follows from (b) since ±||/|| L=°(n) ar e sub- and supersolutions of (3.1). □ 

The limiting problem can be identified through the half relaxed limit method: 

Theorem 3.2. Assume (H/), (H M ), and (H*) for i — 0,1,2,3 hold. Then the 
half-relaxed limit functions 

u(x) — limsup Uk{y) and u{x) = liminf Uk(y) 

are respectively sub- and supersolutions of the Neumann boundary problem in the 
sense of Definition\2.f\ 



In the proof we will need the following result whose proof is given at the end of 
this section. 

Lemma 3.3. Assume (H^) holds for i = 0,1,2, (H M ) holds with c = 1, and let 
5 > and Jn kt r(x) := JLi <r . v( x > z ) dfj. k (z). If i/t^ie dfl as k — > oo, then 

\lu. k ,5{Vk)\ -> °° and -> -n, 

where n = (0, 0, . . . , 0, — 1) is an outward normal vector of dtt. 



Proof of Theorem 3.2 Since the proofs are similar for u and u, we only do the one 
for u. Let S > and <j> € C 2 , and assume that u — <fi has a maximum point x in 
B(x, c v S) D O. Let us first consider the case when x <E fl, i.e. when x^ > 0. By 
modifying the test-function, we may always assume that the maximum is strict. 
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By standard arguments, Uk ~ 4> nas a niaximum point yk in B{x,c rj 5), and when 
k — > +00, 

y k -> x and u k {y k ) u(x). 

Let 5 k = 5 — \x — yie\ and < r < 5k, and note that B{yk, c n r) C B(x, c n 5). Since 
the max of (uk — 4>) in B(y k ,c n r) is attained at yk, we find that 

{I/j. k )r[uk](yk) ■= / uk(yk + v(yk,z))-u k (yk)dfj, k 

J\z\<r 

< / <t>(Vk + viVk, z)) - (t>{yk)dii k = (In h )r[4>](Vk)- 

J\z\<r 



Hence, since Uk is a pointwise solution of (3.1), we find for all < r < 5k, 

u k (yk) - (lp k )r[<f>](yk) - (WKK^fc) < fivk), 

where (7 Mfc ) r [u fc ](a;) := J^ r Uk(x + r](x, z)) - u k (x) d/J, k (z). 

We want to pass to the limit in this equation and consider first the (I^ k ) r -term. 
By the definition of u and (H^), 

limsupu fc (?/fc + rj(y k , z)) < u(x + rj(x, z)) for a.e. z. 

k— >+oo 

Hence, since we integrate away from the singularity of //, we can use Fatou's lemma 
and (H*) and (H^) to show that 

\imsup(I^ k ) r [u k ](y k ) < / u(x + i](x,zj) - u(x) dfj,(z) = I r [u](x). 
k— >oo J\z\>r 

To pass to the limit in the (7^ fc ) r -term, we have to write it as 

<t>{x + v(Vk,z)) - (f>{y k ) - D(f>(y k ) ■ rj(y k , z) dfi k (z) +^ llk . r {y k ) ■ D(/>(y k ), 

\z\<r 



=(I» k )r[<t>](Vk) 

where 7 Mfc ,r(^) := Ji z i <r vi x , z ) d/J, k (z). For \z\ < r, a Taylor expansion then yields 

\HVk +v(Vk,z)) - (f>{y k ) - D(f>(y k ) ■ v(Vk,z)\ < ||-D 2 0||i,«,(B(x,c,r)) \z\ 2 . 

Hence by (H*), (H?), (H*) and (H^), we can use the Dominated Convergence 
Theorem to show that 



(VVMO/fc) -> / (f>(x + rj(x, z)) - <j>{x) - t](x, z)D<j>(x) dfi(z) = I r [(f>](x). 

J\z\<r 

Next, by Lemma [2. 1| 



l^ k AVk)= \ v(Vk, z) d/i#(z) + c r)(y k ,z)dfx k (z), 

where the last integral is understood to be zero if y k ,N > t. Note that since 
Vk,N ~^ %n > 0, the domain of integration of the ^-integral is always bounded 
away from z = when k is big. Along with (HjJ) and (H^), this allows us to pass to 
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the limit in the /i* -integral using the Dominated Convergence Theorem. Similarly, 
we may pass to the limit in the /i#-integral by (H*), (Ff?) and (H ). We find that 



l^ k AVk,N) -> lr{x) : = 



z d[i#(x) + c 



\z\<r 



x N <\z\<r 



r](x, z) d^(x) 



P.V. 



r)(x,z)d(i(x) =■ lr{x), 



\z\<r 



where we used Lemma |2.1| again. Hence we can conclude that 

Hm {Iy, k )r[(t>){yk) = Ir[4>]{Vk) +lr{x) • D(j)(x) = I r [(j)]{x). 

k—>oo 

Since 5k —> 8, we end up with the following limit equation, 

u(x) - I r [<f>](x) - I r [u}(x) <f(x) 

for every < r < S. Using the Dominated Convergence Theorem again, we send 
r — > S and obtain the subsolution condition for (1.1 ) at the point x G fl. 

The second part of the proof is to consider the case of when x £ d£l, i.e the 
case when xm = 0. We first do it in the case c = 1. By adding, subtracting, and 
divinding by terms, we may rewrite the subsolution condition as 

u k {y k ) - {I^ k )sW\{Vk) - {I^ h ) 5 [uk\{yk) ~ f{Vk) li* h AVk) ■ Dcj>{y k ) 



\l^AVk)\ 



\i^ k Ayk)\ 



< 0. 



By Lemma 3.3 \^fj, k) s{yk)\ oo, and since Uk and / are uniformly bounded, 



u k {Vk) 



f{Vk 



(l„ k ) 5 ^k](yk) 

\in k Avk)\ ' \j» k Avk)\ ' \7^ k Avk)\ 

all converge to zero. The same is true for 

(i» k )s[4>](yk) 
h^Avk)\ 

since the integrand of the numerator is controlled by \\D <f>\\<x 



3.3 



satisfies (H*). Using Lemma 

that we may go to the limit in the above inquality to find that 

dc/) 



2 1{ |^| <5} and fj, k 

again, we have l^Av^/b^Av^ -> ~ n , so 
lality to 

(x) < . 



oxn on 

In the case when c = 0, the measure /i = which less singular than ji*. The 
same line of arguments as in the proof for x € f2 (much easier this time) now shows 
that the equation holds at x € <9f2. □ 



Proof of Lemma \3.3[ First note that by Lemma 2.2 with y^ instead of x and /j, k 
instead of (j,, 



ln„AVk)' = / v(Vk,z)' d(i k {z) = 



\z\<6 



\z\<6 



v(Vk,z)' dfi#(z) 



which remains uniformly bounded in k because of (H*) and our assumption on n#. 
Since yk.N xn = 0, we can assume that < yk,N < S, and by Lemma 2.2 



il^ANiVk) 



\z\<8 



y k ,N<\z\<s 



v(yk,z) N d^ k (z) 
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As above, the first integral is uniformly bounded as k — > oo. For the second one, 

(zn - Uk,N) d(j.*{z) > 0, 



we send r — > in Lemma 2.2 (m) to find that 

rj(y k ,z) N dfj£(z) > 



(3.2) 



\z\<6 



■4<s 

yk,N < Z N 



and, since y kt N ~~ ^ 0, we can then use Fatou's lemma to show that 
/ z N dji^z) < liminf / (z N - y kN ) dfi^z) 

J M<.5 k— >oo J M<« 



Applying symmetry of the measure /i* twice, we are lead to 



z N dfi„(z) = - 

z\<S Z 
z N >0 



\z\<6 



\z N \ dpL*(z) 



1 

2N 



N 
■ 

\*\< & i=l 



d^(z), 



so by taking (H^) into account, there is a constant C = C(N) > such that 



^Z N d^{z)> — 

J>0 



\z\<S 



\z\ d/i* (z) 



Hence we have proved that (7^,5)^(1/*;) 

ly boundc( 



{lti k ,&)' {Vk) is uniformly bounded, we see that 



as — > 00, and if we use that 



l7/**,«(»fc)l v bw,*(yk)l ' \in k Avk)\ 



(o,o,..- ,0,1) 



□ 



4. Comparison in non-censored cases 



In this section we prove a comparison result for the non-censored cases, i.e. 
under assumption (H^) for i = 0, . . . , 5. These assumptions covers all the examples 
given in the introduction, except example (a) - the censored case. As a conseqence 
of the comparison result and the results of the previous sections, we also obtain 
well-posedness for (1.1 1. The comparison result is the following: 

Theorem 4.1. Assume (H M ), (H/), and (H^) hold for i = 0,1,2,3,4,5. Let u 
be a bounded use sub solution of (1.1) with data f G Cb(M N ), v be a bounded Isc 
super solution of (1.1) with data g £ C(,(R ) such that f < g in fi. Then u < v on 
Q. 

From this result it follows that the half-relaxed limits in Theorem |3.2| satisfy 
u < u in f2. Since the opposite inequality is always satisfied, this means that u := 
u = u is solution of (1.1 1 according to Definition 2.1 Uniqueness and continuous 

by standard arguments and we have 



4.1 



dependence (on /) follows from Theorem 
the following result: 

Corollary 4.2. Under the assumptions of Theorem ^. 1\ there exits a unique vis 
cosity solution of (1.1 1 depending continuously on f. 



Proof of Theorem \4-l\ We argue by contradiction assuming that M := supq(m — 
v) > 0. We provide the full details only when c = 1. The case c = is far simpler 
since the equation then holds even on the boundary. 
To get a contradiction, we first introduce the function 

*r(x) := u(x) - v(x) - i/j r (x,x) , 
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where ij)R is a localisation term which makes the max attained: 

1> R {x, y) = rf>(\x N \/R) + 4>(\y N \/R) + *P(\x'\/R) + tl>(W\/R) , 
with ip a smooth function such that 

!0 for0<s<l/2, 
increasing for 1/2 < s < 1 , 

2(|M| 00 + H 00 + 1) forsM. 

Of course, Mr '■= max'&R — > M as R — > oo so that for R big, max^R > 0. Then 
there are two cases for such R > 0: 

(a) either there exists a maximum point xr of ^r which is not located on the 
boundary. In this case the proof is quite classical: we use the doubling of variables 
plus a localisation term around xr by considering the max of 

I IjJ y ' I 2 I ^J^y U N I 

u(x) - v(y) ^ 2 iPr( x i v) - °\ x - xr\ 4 > 

e e N 

where a > is small. The localization term a\x — xr\ a is chosen so that xr is 
the unique maximum point of x i— > ^r(x) — cr|x — xr\ 4 , and by choosing a small 
enough, the contribution of this function in the integral term is smalQ 

Hence the maximum points (5, y) of the above test-function converges nccessarely 
to {xr,xr) as e',£n and they are also bounded away from the boundary if 
e', En are small enough. This property implies that we can use directly the equation 
and obtain max ^ r < 0, which is a contradiction. 

(b) or any maximum point xr is located on the boundary. In this case we use the 
doubling of variables plus some extra term to push the points inside (see below) 1 : 

®e>,e N ,u,R( X >y) := u ( x )~ v (y)-— W~ ~ 2^" V^O, y)+d„(x N )+d„{y N ) . 

In this case we can assume without loss of generality that the maximum points x, y 
are always such that < xn, JJn < 1, whatever ejv, e', v > are. 

Note that in both cases we take two distinct real parameters e',£n > in order 
to take into account the special contraction property in the iV-th direction, see 
(H^). Now, since case (a) is rather standard, we only provide a proof of case (b) 
which is more involved. 

The term d^ plays the role of a distance to the boundary of the domain; such 
term is usual in classical Neumann proofs in order to prevent the maximum points 
to be on the boundary. More precisely, for v > 0, we take d„(-) = ^d(-) where d is 
a smooth function such that 



d(s) = { 



s for0<s<l/2, 
increasing for 1/2 < s < 1 , 
1 for s > 1 . 



Let us note that if < v < 1 and R ^> 1 are fixed, then $ := <& e ',E N ,v,R < 
for \x\, \y\ large enough, while, by choosing x = y in a suitable way by taking into 
account the fact that M > 0, we have $ e < ejv v r{x,x) > for v small; hence the 



^In the viscosity inequalities, the various penalization terms are only integrated near the origin, 
e.g. in {\z\ < S}. Therefore we do not need to worry about the integrability at infinity of \x\ 2 and 
|i| 4 with respect to the measure fi. 
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maximum of $ is attained at some point (x,y) £ fi 2 , that we still denote by (x,y) 
for simplicity. 

After proving that the points x, y are inside f2, we are going to first let ejv — > 0, 
then e' — > 0, then v — > and finally R oo. Because of this use of parameters, we 
have 

M SN>s , tI/jR := max$ £jv , e>!fl -> M > 0. 
In particular, this implies that 

|x' - y'l 2 

x w - y N = 0(e N ) , x' - y' = O(e') , — — — = o £N)£ /(l) , 

where the 0(eAr),0(e') are uniform with respect to all the parameters, and the 
o EJViE /(l) means precisely that after passing to the limit as Sn 0, we are left with 
a quantity which is an o e /(l). Also note that 

u(x) - v(y) = M + o £N ^. v ^ R (l) , 

where the order of the parameters is important as explained above. 

Step 1 — Pushing the points inside. 
We denote by 

I \ \ X '-V'\ 2 , \%N-yN\ 2 . , , \ A ( \ A r \ 

e e N 

where we have dropped the parameters for the sake of simplicity of notations. 

In this step, we prove that the i^-viscosity inequalities hold for u and v. Accord- 
ing to Definition [O] this is clearly the case if c = since these viscosity inequalities 
hold even if the maximum or minimum points are on the boundary. 

In the c = 1 case, let us assume that the maximum point (x,y) is such that 
xyv = 0, then a; is a (globa l) m aximum point of the function z i— > u(z) — v(y) — <p(z, y) 



and, thanks to Definition 2.1 we should have —-§^{x,y) < 0. But, recalling that 



is zero in a neighborhood of the boundary, we have 

dip 2(x N - y N ) dip R d 2y 



., 0,2/) = — ~ 2 - ^^(x,y) + - r d^(0) = -~ +v > , 
oxn £ n oxn ds e z N 

which is a contradiction. Therefore xn cannot be zero and a similar argument 
shows that yN > as well, hence both x and y are inside Q. 

Step 2 — Writing the viscosity inequalities and sending 5 to zero. 
We introduce a (small) fixed parameter < 5 < 1 which is the parameter appearing 
in Definition [2?T] in order to give sense to different terms in the equation. We write 
the definition of viscosity sub and super solutions, using the test-function in the 
ball B$ for a S < p := min^N, y^, 1), and the functions u and v outside this ball. 
Since u is a viscosity subsolution and the function u(-) — v(y) — <p(-,y) reaches a 
maximum at x, then we have the viscosity subsolution condition that we write as 
follows, thanks to Lemma 12. II 



u{x) - I [<p(x + r)(x, z), y) - ip(x, y) - D x <p(x, y) ■ rj(x, z)] dfi(z) 

J\z\<& 

-P.V. / D x <p(x,y)r)(x,z)dfj,(z) - [u(x + rj(x,z)) - u(x)]dfj,(z) < f(x) . 

J\z\<8 J\z\>6 
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For simplicity of notations, we leave out the P.V. notation since the integral can 
be expressed as converging integrals and we use the notation P(x, z) := x + rj(x, z). 
Next we use Lemma |2.2|-(z) and the similar super solution condition on v to get 



[<p(P(x, z),y)- ip(x, y) - D x tp(x, y) ■ t](x, z)] dfi{z) 

\z\<S 

D x tp(x,y) ■ / r]{x, z) dfj,#(z) 

J\z\<8 

[ip(x, P{y, z)) - ip(x, y) + D y tp(x, y) ■ r){y, z)] d/j,(z) 

\z\<S 



+ D y cp(x,y)- rj(y,z)dn#(z) 

J\z\<S 



[u(P(x, z)) - v(P(y, z)) - u(x) + v(y)} d/J,(z) 

l\z\>S 

+ u(x) - v(y) < f(x) - f(y) . 
In order to pass to the limit as 5 — > to get rid of the test-function tp, we use 



Lemma 2.1 for all the terms which are smooth functions: the integrals over B(0, 5) 
all vanish as 8 — > and we are left with limit of the integral over {\z\ > 8}. To 
this end, we split this integral into two integrals, one over {|z| > 1} (which is 
independent of 8 of course) and the other over {8 < \z\ < 1} that we have to deal 
with. 

Using the definition of the maximum point for $, we have that for any z: 

u(P(x, z)) - v(P(y, z)) - p(P(x, z), P(y, z)) < u(x) - v(y) - ip(x, y) . 
Hence, it follows that 

u(P[x, z)) - v(P(y, z)) - (u(x) - v(y)) 

< \P(x,z) N ~P(y,z) N \2 \x N -y N \ 2 | \P(x, z)' ~ P(y, z)f \x>-y>f 

+ i>B.(P{x, z), P(y, z)) - ip R (x, y) 

- d u (P(x, z) N ) + d v (x N ) - d v (P(y, z) N ) + d ly (y A r) , 

and we put this inequality into the integral over {<5 < \z\ < 1} which gives rise to 
several terms denoted by (with obvious notation): 

{u(P(x, z)) - v(P(y, z)) - u(x) + v(y)} dfx(z) < T s Sn + T s e , + T* R + . 

8<\z\<l 

As for the £Ar-terms, we get rid of them by (Hjj) which implies T^ N < 0. Then for 
the e'-terms we write 

T s = f , \p(x, z y-p(y, z y\* W-y'\ 2 \ dfi 

£ Js<\z\<l ^ £ ' 2 £ ' 2 ' 

<\ I \v(x,z)' - f](y,z)'\ 2 dn(z) 

£ J<5<|z|<l 

(x' - y') ■ (t](x, z)' - T](y, z)') dn{z) . 

8<\z\<l 
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For the first term of , we use the domination of the integrand by c\z\ 2 to pass 
to the limit as 5 — > 0. For the second one, we use Lemma 2.2 -(Hi) which allows to 
wipe out the symmetric /i* -contribution, so that we get in the limit 



lim sup T*, < —r 
5-»o £ 



\r)(x,z)' -r}(y,z)'\ dfi(z) 



0<|z|<l 



2 

7^2 



{%' - y') ■ z)' - r)(y, z)') dfj,#(z) . 

>0<\z\<l 

We concentrate now on the penalisation terms which are given by integrals of 
smooth functions. Note first that we are in the case when < xn,Un < 1, so that 
the iI>( x n/R) and ipiljN / -ff)-terms vanish (we assumed that R^> 1). Hence, using 
Lemma [271] we get as 5 — > the following two contributions: 



lim Tl 

5^0 



5 



hK](x)- -d u (x N )P.V. 
ds 



r)(x,z) N diJ,(z) + (...)(y), 



0<|z|<l 



limT* s < - h$ R ](x) - Dj> R {x)' ■ P.V. 



0<|z|<l 



■q{x,z) N dfi(z) + (... X^- 



where ip R (x) := ip(\x'\/R) and the (. . . )(y) notation stands for the same t erms but 
calculated at y instead of x. Now, note that 4-d(x^) > and use Lemma 2.2-(m) 
(with r = xn > 0). This gives that in the principal value for T^ ul the /^-term 
which is multiplied by {—v) has a nonpositive contribution. So we find that 



lim T? < -vl^dUx) - v — d{x N ) 
<5->o " ds 



0<|z|<l 



7/(35, z) N d^ # {z) + (...)(y) = o„(l) . 



As for the T^-term, this time we use Lemma 2.2 (zi), which implies that the 
symmetric /i^-part of the principal value vanishes: 



lta^ = -7 1 W(i) 



1 

R 



<C(p) 



1 



-(\x'\/R) 
\i>\\c 



'0<|z|<l 
Ofl(l). 



Tl(x,z)' dn # (z) + (. . .)(y) 



ds 

R2 l\f lie R < 

Thus, the parameters e',sjsr,u,R > are still fixed for the moment and after 
sending <5 — > we have obtained: 

u(x) - v(y) < f(x) - f(y) + o„(l) + o R (l) 

1 



\r)(x,z)' -r){y,z)'\ dfi(z) 



\z\<l 



(x' - y') ■ (<q(x, z)' - r)(y, z)') dfi # (z) 

\z\<l 

+ / {u(P(x,z))-v(P(y,z))-u(x) + v(y)}dfi(z) 

J\z\>l 

= f(x) - f(y) + o v {l) + o R (l) + Inti + Int 2 + Int 3 . 



Step 3 — Sending the parameters to their limits 

We let first ~~ ► 0, the other parameters remaining fixed for the moment and we 
recall that \x N — y N \ = 0{en). Moreover, as long as R is fixed, the points x,y 
remain in a compact subset of fi; therefore we can assume without loss of generality 
that x, y are converging to points (still denoted by x,y) such that x^ — y^. 
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Combining (H?) and (H^), we obtain in particular that 

lim \r)(x, z)' — rj(y, z) \ < C\z\\x — y \ for fi-a.e. z. 

Then, (H*) and the integrability condition on justify that we can use dominated 
convergence in Int2. The argument is similar for Inti, using the domination 

\ V (x,zY-r 1 (y,zy\ 2 <(2c v ) 2 \z\ 2 . 

So we find that lim EJV _>o Int2 = while 



lim Inti < C ^— [ \z\ 2 dfi(z) = o £ >(l) 

e J\z\<i 



The Ofl(l) and o v (l) terms are uniform with respect to the other parameters, so 
there is no problem to send £jy, e' — > here. Next, since \x — y\ — > as s^, e' — > 
here, we may assume that x, y — > x € f2 by considering a subsequence is necessary. 
By continuity of /, it then follows that (f(x) — f(y)j —> 0. 

We then pass also to the limit as v — > and get: 

(4.1) u(x) — v(x) < limsup limsup limsup [Int 3 ] + <3r(1) . 

v-yQ e'-H1 s N ->0 

Passage to the limit in the Int 3 term is possible because the domain of integration 
does not meet the singularity of the integral: we need only use the u.s.c. and 
l.s.c. properties of u and v, together with Fatou's Lemma (because the integrand 
is bounded and fx is finite on {\z\ > 1}). After passing to the limit in En,£' and v, 
we have by definition 

lim lim lim (u(x) — v(y)) = M + <3r(1) 

v— >0 e' — >0 ejv— >0 



limsup limsup limsup [Inta] 

u->0 e'->0 £n->-0 



so that 



< / {u(P(x,z))-v(P(x,z))- (M + o R (l))\ d(j,{z) . 

J\z\>l 

Now since and u(P(x, z)) — v(P(x, z)) < sup^(u — v) = M, 

limsup limsup limsup [lnt^\ < / °fl(l) d/j,(z) = or(1). 
When R oo in (4.1 1, we get M < and the proof is complete. □ 

5. Comparison in the censored case I. 

In this section we give comparison and well-posedness results for the initial value 
problem (1.1 1 in the censored case (under assumption (Hf ; )) when the measure // 
is not too singular: 

(H£J The measure /i is a nonnegative Radon measure satisfying 



(i) / \z\ A 1 d[i < oo and (ii) / d[i = for any a < 

JR N J{z N =a} 

In addition, we assume the existence of a "blow-up supersolution" 
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(U) There exists Rq > such that, for any R > Ro, there is a positive function 
U R £ C 2 (Q) such that 

-I[U R ]{x) > -K R in {x : < x N < R}, 
for some Kr > 0, and 

U R (x) > T — r in il, 

ur(x n ) 

for some function lor which is nonnegative, continuous, strictly increasing 
in a neighbourhood of 0, and satisfies w(0) = 0. 

Remark 5.1. See Appendix [A] for a discussion of this assumption. E.g. in Remark 



A.l we prove that (U) holds if 

M 

ij, = p, + 2J Ci§ x i 

i=l 

where a £ K, S x i are delta measures supported at {x 1 } for x l N > 0, and 

^ = J^L where «e(0,l), 0< 5 eL°°(H), limg(z)=g(0)>0. 

This class of measures include the Levy measures of the stable, tempered stable, 
and self-decomposable Levy processes. Much more general examples are presented 
in Appendix [A"| 

Theorem 5.1. Assume (HJJ, (H/) ; (H^j) and (U) hold. Let u be a bounded use 
subsolution of (1.1) and v be a bounded Isc supersolution of ( |1.1[ ). Then u < v in 

n. 



As in the previous section, we immediatly get a well-posedness result for (1.1) 



by Theorems 5.1 and 3.2 



Corollary 5.2. Under the assumptions of Theorem \5.1\ there exists a unique vis 
cosity solution of (1.1) depending continuously on f. 



Proof of Theorem \5.1\ We argue by contradiction assuming that M :— supjj (u — 
v) > 0. Take R > Rq and < k <C 1. Using < e <C 1, we double the variables by 
introducing the quantities 

I 1 2 

4>(x, y) = ]X ~ Vl + k[U r (x) + U R {y)] + 4> R {x) + ^ R (y), 



$(x, y) = u(x) - v(y) - cf>{x, y), 

.,11 _i_ IL, II 

R 



where Ur is given by (U) and iPr(x) = 2(||u||oo + ll w l|c»)'0(^) f° r an increasing 



function t/j(s) £ C°°(0, oo) which is in (0, \) and 1 in (1, oo). 

For any R, k and e, the function $ achieves its maximum at (x, y) = {xr^ K}E , 
and, by the definition of Ur and ipR, we have 

(5.1) XN,m>do = u)n 1 ( J and \x\,\y\ < R. 

\2(\\u\\ 00 + WvWoo)/ 

These estimates will hold in most of the proof since we are going to keep R and n 
fixed untill the end, sending e —¥ first. A standard argument also shows that 
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By the estimates on x, y and extracting a subsequence if necessary, we can assume 
without loss of generality that x,y —¥ X, u{x) — > u(X), and v(y) —> v(X) where X 
is a maximum point of Q(x, x) = u(x) — v(x) — 4>{x 1 x). Finally, when we first send 
k — > and then R — > +00, we have 

u(X) - v(X) ->■ M and kUr(X) + ip R [X) -> . 

Now we write down the viscosity inequalities. Since u — (f>(-,y) has a global 
maximum at x and v — (— 4>(x, •)) nas a global minimum at y, we have that 

u(x) - I 5 [u}(x) - l s [4>(.,y)](x) < f(x), 

v{y)~I s [v]{y)~h[^{x,-)](y)>f{y). 

With this in mind we see that 

M + o(l) = u(x) - v(y) - cf)(x, y) 

(5.2) < I s [u]{x) - I s [v](y) + I 5 [4>(; y)](x) - h[<j>(x, -)]{y) + f(x) - f(y). 

In this inequality, we aim at first sending 6 — > in order to get rid of the e-depending 
7,5 [</>] -terms. In fact Is[f] — > as S — > by the Dominated Convergence Theorem 
since |7y(a;,z)| < c n \z\, and hence for any C^-function tp, 



\cp(x + 7](x, z)) - ip{x)\ l\ A<5 d^{z) < Cri \\Dip\\ Lx(B ) I l\ z \ <6 \z\ dn(z). 

Next we consider the /''-terms. We restrict ourselves to a subsequence such that 
%n > yN (if %n < Hn the argument is similar). Then 

I s [u](x) - I s [v}(y) = / [u(x + z)- u(x)] l lz>s dn(z) 

J -x N <z N <-y N 

[u(x + z)- v(y + z)- (u(x) - v(y))} l\ z i >s diJ,(z) 

-Vn<Zn 
h+h- 



For Ji, we have 



\h\ < 2||m||oo / l { - SN<ZN< ^y N} (z) dn(z) 

J\z\>& 



<\z\>8 

Keeping n and R fixed and recalling (5.1 ), we see that this integral is independent of 
5 as soon as S < S . Furthermore, because of (HJJ (ii), the Dominated Convergence 
Theorem implies that 

1 1 ->• as e ->• 

since \x — y\ — > as e — > 0. 

For J2, we use the maximum point property for x, y, 

(u(x + z)- v(y + z)) - (u(x) - v(y)) < <p(x + z,y + z) - (j>(x, y) , 

which after cancellation of the e-terms leads to 

h < K (l 5 [U R ](x)+I s [U R ](y)) + (l 5 bP R }(x) +I S [^ R }(y)) ■ 

Recalling again ( |5.1[ ) and using the regularity of U R and <p, we can send 8 — > and 
obtain 

hmsup I 2 < k(i[U r ]{x) + I[U R }{y)) + (l[il> R ](x) + I[ip R ](yj) , 
where each term on the right-hand side have a sense. 
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Consider equation (5.2) again. Using all the previous estimates, we can send 



<5—^0 first and obtain using (U) for the W^-terms that 

M + o(l) < 2K r k + (I[iPr](x) + I[ip R ](y)) + (f(x) - f(y)) . 

In this inequality, we can first send e — > 0, keeping R and k fixed. Then f(x) — 
f(y) — > as e — > since / is uniformly continuous in Br, and we find that 

M + o(l) < 2K r k + 2I[i> R ]{X) . 

We conclude by first sending k — > and then R — > +oo. □ 

6. Comparison results in the censored case II. 
In this section we give comparison and well-posedness results for the initial value 



problem (1.1) in the censored case (under assumption (H^)) when the measure /x 
is very singular 

(R'P Hypothesis (H M ) holds with 

/i*(rfz) = j m+a , [ (1 A |z|^) (i#(dz) < oo , f /i^(dz) = for any a < , 



{z N =a} 



for a E (1, 2) and j3 := a - 1. 

This assumption is much more restrictive than (H M ), and the results of this sec- 
tion are not completely satisfactory. We had lot of difficulties to obtain comparison 
results because on one hand, it is not possible to get rid of the boundary and the 
boundary condition in such a general way as we did in the less singular case I. On 
the other hand a lot of technical difficulties come from the the way the x-depending 
domain of integration in / interferes with the singularity of the measure and the 
boundary. 

Our first result is the following 

Theorem 6.1. Assume (TL f ), (H^), and (IT7) hold. 

(a) Let u and v be respectively a bounded use subsolution and a bounded Isc super- 
solution of 

(6.1) w(x) — I[w](x) = f(x) in Q , 

and let us also denote by u and v respectively their use or Isc extensions to ?Q // 
there exists C > and ft > ft such that 

(6.2) u(x' ,xn) > u(x' ,0) — Cx^ and v(x' ,xn) < v(x' ,0) + Cx 1 ^ 

then u and v are respectively a bounded use subsolution and a bounded Isc superso- 
lution of 

(b) If u and v are respectively a bounded use subsolution and a bounded Isc super- 
solution of ( |1.1[ ) satisfying (6.2), the 



len 

u < v in fl. 



In particular, there exists at most one solution of ( |1.1[ ) in C 0, ^(J7) for (5 > (3. 



2 For any x' £ H N 1 ,u(x',0):= limsup u(y' , i/jv) and v(x' , 0) := liminf f(y',3/jv) 

(a',!/jv)->(a:',0) (»',!«v)-Kx',0) 
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Several comments have to be made on the different statements in Theorem iG.f I 
Part (a) means that, for sub and supersolutions having a suitable regularity at the 
boundary, the Neumann boundary condition is already encoded in the equation 
inside. This might be expected from the proof of Theorem |3.2| or from the intuition 
coming from the censored process. But the result is not true in general since we 



need anyway (6.2) to prove it. 

Unfortunately part (b) does not provide the full comparison result for semi- 
continuous solutions, and we do not know if this result is optimal or not. Of 



course, in view of Theorem 6.1 (b), it is clear that we need a companion existence 



result providing the existence of solutions satisfying (6.2 1 or be long ing to C 0, ^(f2) 



for j3 > (3. We address this question after the proof of Theorem 



6.1 



Proof. We prove (a) only in the subsolution case since the supersolution case is 
analogous. Let <j) be a smooth function which is bounded and has bounded first and 
second-order derivatives and assume that u— <j> has a maximum point (x' , 0) € dtt in 
B((x' , 0), c v 5) n ft. We may assume that the maximum is strict and global without 
any loss of generality. 

We set 9{t) = t 13 A 1 for t > and, for < n <C 1, we consider the function 
u(x) — 4>{x) + k6(xn). By standard arguments, using the properties of 0, this 
function achieves a global maximum at a point nearby (a/,0), and we claim that 
this point cannot be on dtt = {x : xn = 0}. Indeed, otherwise it would have to be 
(x' , 0), the strict global maximum point of u — <fi on dVt. But then by ( |6.2[ ), 

u(x', 0) - 4>{x', 0) > u{x) - 4>(x) + k9(x n ) > u(x', 0) - <f>(x', 0) - 2Cx% + k9(x n ), 

and we have a contradiction since (3 > f3 and hence —2Cx N + k&(xn) > for xn 
small enough. 

Therefore the function x n- u(x) — 4>{x) + k8(xn) has a maximum point at x K 
with (x k )n > 0. We may write the viscosity inequality at 

u(x K ) - I s [<p](x K ) - j{x K ) ■ Dcf>(x K ) + kI s [9](x k ) - I s [u](x K ) < f(x K ), 

for (say) < 5 < 1, where j(x K ) = P.V. f^ <g r](x K , z)fi(dz). 

We first consider the term kI$[9](x k ). On one hand, the /i^-part is O(n) since 
9 is in C '^ and (H^) holds. On the other hand, the singular part (the //* part) is 
nothing but 

f dz 
kF Y J k ,< s ( x n + zn)-9(x n ) , 

where we have dropped the subscript k to simplify the notation. Since 5 < 1 and 
ijv — > as K — > 0, we may assume that < xn + zjv < 1 for \z\ < 5, and hence 
that the principal value reduces to 

KYN.j \x N + z N f-\x N f-^. 

By the computations of Lemma |B.1| in the Appendix, 



— — rl7 

-P.V. / \ XN + ZN f-\ XN f—— = Vl 
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for xn > 0. Writing 

«P.V. / (•••; '. I'.v. / (•••) -k f (■■■), 

we conclude that for fixed 5, 

f dz 
kP.Y. 6(x n + z n ) - 6(x n ) -—- = 0(k). 

J *N+*N>° |Z ' 

Finally, the u, Is, and I s terms are uniformly bounded in k while j(x K ) oo 
since (x k )n — ► 0. We divide the above inequality by |7(x K )| and send k 0. As in 
the proof of Theorem 3.2 - the second part, when x € dft and c = 1 the result is 
that all terms vanish except the 7-term and we are left with the boundary condition 

Now we prove part (b). By linearity of the problem and part (a), the function 
w = u — v is a subsolution of with / = 0, and we are done if we can prove 
that w < 0. To prove this, we consider the function 

XrA x ) ■= tp(\x N \/R) + i>{\x'\/R) - vd(x N ) , 
where ip and d are defined as in the proof of Theorem |4.1[ replacing, in the case 
of ip, 2(\\u\\oa + |H|co + 1) by 2||iy|| 00 + 1. The function xr.v is smooth and easy 
computations show that xr,v is a supersolution of with an / > zu(R, v) where 
zu(R, v) converges uniformly to as R — > 00 and v — > 0. At the boundary <9f2, 

_^ =0 + ,.l>0. 

OXn 

Because of the behavior of Xr,v &t infinity, the function w — Xr,v achieves its 
maximum at some point x, and because of the behaviour of XR,v a t the boundary, 
xn > 0. Writing the viscosity subsolution inequality then yields that 

w(x) - XrA x ) ^ ~XrA x ) + i[xrA{ x ) + 1& \ u ~ XR,v]( x ) ^ -™(R, v) + 0, 
where we have used that I s {ip](x) < at any maximum point x of ip. Hence, for 
any y G fi, 

w(v) - xrAv) ^ -™(R, v), 

and part (b) follows from sending R — > 00 and then v — » 0. □ 

Now we turn to the existence of Holder continuous solutions and we begin with 
a result in 1-d. 

Theorem 6.2. Assume N = 1 and that (H/), (H®), and (H^) hold. 

(a) Any bounded, uniformly continuous solution of (1.1) is in C '^(O) for some 
f3>p. 

(b) There exists a solution of (1.1) in C 0, ^(f2) for some (3 > (3. 
Proof, (a) To prove the Holder reglarity we consider 

(6.3) M= sup (u(x) - u(y) - C\x - yf) , 

[0,+co)x[0 : +oo) 

and argue by contradiction assuming that M > 0. The aim is to show that this is 
impossible for C > large enough. A rigorous proof would consists in introducing 
localization terms like the -0-terms in the proof of Theorem |4.1| and d^-terms in 
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order to take care of the Neumann boundary condition. We drop these terms for 
the sake of simplicity in order to emphasize the main idea and not loose the reader 
in technicalities. 

Therefore we assume that the above supremum is achieved at (x, y) with x, y > 0. 
Since M > we have x ^ y, and we assume below x < y. The other case 
can be treated analogously. To simplify the notation, we introduce the function 
4>(z) :— C\x — y + z\P . Note that this function is concave in the intervals (-co, y — x) 
and (y — x, +oo), and that it is smooth in (—6, S) for 5 < y — x so that it can be 
used as a test function. By the maximum point property for (x,y), 

u(x + z\) - u(y + z 2 ) - C\x — y + (zi — z 2 )\ 13 < u(x) - u(y) - C\x - yf , 

for z\ > —x and z 2 > —y, and hence 

(6.4) u(x + z) ~ u(y + z) — [u(x) — u(y)] < for z > — x (> —y), 

(6.5) u(x + z) — uix) < [4>{z) — 0(0)] for z > —x, 

(6.6) u(y + z) — u{y) > — [<j>{— z) — 0(0)] for z>—y. 

Using the definition of viscosity solution and the symmetry of the measure /x* , 
for 6, 6' > small enough, we have the inequalities 

-(Is[^+I 5 [u})(x) + u(x) <f(x) and - (!*#] + I S '[u])(y) + u(y) > f(y), 
which reduce here to 

— 6 p5 

(u(x + z)- u(x))dfx(z) - J [<j){z) - (f>(0) - 4>'{Q)z]dn{z) 

5 r+oo 

(jj (0)zdfi(z) — / (u(x + z) — u(x))dfj,(z) + u(x) < f(x) , 

-S J 6 

(u(y + z)- u(y))dn(z) + / [</>(- z) - 0(0) + (0)z]dn{z) 
-y J-S> 

8' p-\-oo 

4>'{Q)zd^{z) - / (u(y + z)- u(y))d»(z) + u(y) > f(y) . 

-8' J 8' 

In the proof below we will subtract these inequalities and the main difficulty of 
the proof will come from the term 

J ■= - / (u(y + z)- u(y))d(j,(z) 
J-y 



which is not a difference of terms from (6.7) and (6.8). Indeed the domain of 
integration z £ (— y, — x) appears in inequality (6.8) but not in (6.7). Because of 
the singularity of fx, if x is close to it is not obvious how to get an estimate for 
J which is independent of C, or how to control this "bad term" by a good term. 
Therefore we have problems with this term if x — > when C — > +oo. For the 
/*t#-part of J there is no problem, we can use (6.6) to see that 

c p— x r~ x 

{u{y + z)-u(y))dn # {z)< [cf>{- z) - </>{0)]dfi # (z) < C \zfd^(z) , 

J-y J-y 

and we will see later that this term can be controlled since \z\@ is /z#-integrable. 
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First case - We first consider the case when x < y — x, or equivalently, 2x < y. 



In this case J > and can be dropped from inequality (6.8). To see this we note 
that for — y < z < —x, 

2.t — y < x — y — z < x 



with x < y — x and 2x — y = — (y — x) + x > — (y — x), and hence by (6.6) 
(6.9) u(y + z)- u{y) > -[(/>(- z) - 0(0)] = \x - yf - \x - y - zf > . 

In this first case, we choose S — x and 5' = y — x and subtract the viscosity 



inequalities (6.7) and (6.8). After some computations using (6.5), (6.6), and (6.9), 



and dropping the J term, we are lead to the inequality 



[<f>(z) + <p(-z) - 2(p(0)]dfj,(z) 



—X 



((u(x + z)- u(y + z)) - (u(x) - u(y)))dn(z) + u(x) - u(y) < f(x) - f(y) 



Some easy computations then shows that the first integral equals 



-C(y-x) 



dz 



\l + g f + \l-zf -2]j^ + 0(C) , 



where the 0(C)-term comes from the \i# part of the measure since the integrand 
can be estimated by 2|z|' 3 which is integrable on, say, (—1,1). The second integral 
is nonpositive by (6.4) and can be dropped because of the "— " in front. 



Finally, since / is bounded and u(x) — u{y) > (by assumption), we obtain 
.10) -C(y-xf~ a f> ■ ' ' ■ - '' : 



[|1 + zf + |1 - zf 2]^ < 2||/|U + 0(C) . 



In order to conclude, we use that M — u(x)~u(y)~-C\x~y\ l3 > (by assumption) 
and j3 < 1 < a to find that 

\x-y\ < f 2 l^y //3 and C(y - xf~ a > K & , 



C J 

where Q:=l + (a- P)^ 1 > 1 and K = (2\\u\ 



J 1 _^[\l + zf + \l-zf-2]- 



dz 



> - 



\l + z\ 



Then we note that 
'' + \l-zf-2] v ^->Q 

I I J U 1+Q 



since z i— > |1 + z\@ is strictly concave on (—1,1). From inequality (6.10) we then 
find that 

KC C < 2H/IU + 0(C) , 
which cannot hold for C large enough and we have a contradiction in the first case. 



Second case - When x > y — x, or equivalently, 2x > y. In this case we choose 
5 = 8' = y — x, subtract viscosity inequalities (6.7) and (6.8), and use (6.6) to see 
that 



[<t>{-z) - cj>(Q)W(z) - I [<j>{z) + <j)(-z) - 2</>(0)}dn(z) 
a J-(y-x) 
r+oo 

I ((u(x + z)- u(y + z)) - (u(x) - u(y)))dfi(z) + u(x) - u(y) < f(x) - f(y) 

y-x 
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Arguing as in the first case, we can drop all it-terms and are lead to an inequality 
of the form 



-C(y - xf- a (B(a) + G) < 2||/|U + 0(C) 



where 



B(a)= (\l + z\ 

J-a-l 
1 

P 



dz 



G= (|i + z|P + |i-^_2) 



dz 



-l 



\N+a ■ 



with a = x/(y — x) > 1. A technical computation (Corollary B.3 in the appendix) 
then shows that B(a) + G < —k < for some (3 > f3 and we can conclude the 
argument as in the first case. The proof of (a) is complete. 

Note the important estimate, valid in both cases: there exist k\,k 2 > such 
that 

(6.11) I[u](x) - I[u](y) < -k x C\x - yf- a + fc 2 (l + C) , 

where the 1 comes from the localization terms. This formal estimate should be 
interpreted in the viscosity sense and with the above choice(s) of test function and 
parameters 5 and 5', cf. e.g. (6.10). 

(b) To show the existence of solutions with a suitable regularity property, we follow 
the so-called "Sirtaki method in 4 steps" . We just give a formal sketch the proof 
which is an easy adaptation of the above arguments. 

We start by building a suitable approximate problem. We approximate the Levy 
measure fi by bounded measures fj, n = /il| z |>i/n for n > 1 and denote the associated 
nonlocal term by /„. Then we introduce a truncation of the nonlocal term and add 
an additional viscosity term. The results is the approximate equation 

-eu xx - T R (I n [u]) +u = f in , 

where Tr(s) := max(min(s, R), —R), R, e > . 

1. For fixed e,n,R, since the T^-term and the measure /i n are bounded, this 
equation can easily be solved by classical viscosity solutions' methods (Perron's 
method and comparison result) . This provides us with a continuous solution which 
is bounded and we even have ||u||oo < ll/lloo- 

Moreover, in order to obtain the C '^ -regularity and C 0, ^-bounds, we consider 
(6.3) and follow the arguments in the first part of this proof. After subtracting the 
viscosity sub- and supersolution inequalities, we formally obtain 

-e [u xx (x) - u xx (y)} - [T R (I n [u])(x) - T R (I n [u])(y)] 

(6.12) +u(x) - u(y) < f(x) - f(y) . 



For the second-derivatives, we have an analogue estimate to (6.11), namely there 
exists k[ , k' 2 > such that 

(6.13) u xx {x) - u xx (y) < -k[C\x - y\^ 2 + k' 2 . 

Note that to give meaning to this formal estimate, we must consider instead of u xx 
the sub- and super jets of the theorem of sums, cf. e.g. Now consider (6.12) 
with fixed R, e > 0. Since the Tfj-terms are bounded, we can rewrite it as 

{y)\ < 2i? + 2(| \u\ loo + 11/Hoo) , 
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and use (6.13) to find that the inequality cannot hold for C large enough. This 
implies that the solution {u n ' R ' e } is at least C 0, ° by the arguments of the regularity 
proof above. 

2. The above argument also shows that, for fixed e, the C^-bounds for the {u n ' Re } 
are uniform in n since they depend only on R through the Tfj-term. This allows 
us to pass to the limit n — > +oo and get a solution u R ' e :— lim„_ i . +00 u n ' R ' e of 
the limiting equation enjoying the same C 0,,3 -bound. This solution satisfies the 
truncated viscous equation with fi n replaced by the singular measure \x. 

3. Next, we repeat the proof of the C^-bound for the truncated viscous equation : 



Estimate (6.11) together with the fact that T R is an increasing and a 1-Lipschitz 

T R (I[u](y))<k 2 . 



continuous function, implies that 

Tr(I[u](x)) 



at least for C big enough. Rewriting the analogue of (6.12) as 

-e[u xx {x) - u xx {y)] < [T R (I[u])(x) - T R (I[u])(y)] + 2(||«|| 0o - 

this new estimates on the difference of the truncated terms shows that the C° : ^- 
bound which is obtained in Step 1, is independent of R and we can let R — > +oo. 
The result is that the limit u e := limfl_ J . 00 u R ' £ is a C 0, ^-solution of the non- 
truncated viscous equation 



-I[u] 



f in n 



4. Finally we come back again to the proof of the C 0, ^-bound but, this time, the 



main role is played by the non-local term via estimate (6.11). Indeed we rewrite 



the analogue of (6.12) as 



- [I[u]{x) - I[u](y)] < e \u xx {x) - u xx (y)] + 2(\\u\ 



and remark that, since the it^-terms satisfy (6.13), the e-term in ( |6.12 ) can be 
estimated by ek' 2 . Using (6.11), we obtain again a contradiction for large enough 
C. The argument is the same as in Step 3 with the roles of the local and nonlocal 
terms exchanged. This also explains the terminology "Sirtaki's method", since 
Sirtaki is a danse where we exchange the roles of the two feet as we exchange here 
the role of the eu xx and I[u] terms. To conclude the argument, we have found that 
the C°'°-bound is independent of e, and we pass to the limit as e — > 0. We get a 
solution u of the original problem belonging to C '@ . Since this solution is unique, 
it is the solution we are looking for. □ 



Now we turn to the case when N > 2. 
retrictive assumptions on /. 



Unfortunately we require far more 



Theorem 6.3. Assume N > 2, that (R f ), (H^) ; and (H'^) hold, and that /(..., x N ) 
is in tU 2 ' 00 (IR Ar_1 ) for any xjy > with uniformly bounded W 2 '°° -norms. 



(a) Any bounded, uniformly continuous solution of (1.1) is in C°'^(fl) for some 



(b) There exists a solution of (1.1) in C°'^(p,) for some (3 > (3 



Proof. We are not going to provide the full proof since it is rather long and tedious 
and is mostly based on two ingredients which we have already seen. But we remark 
that an easy consequence of the the comparison result and linearity of the problem, 
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is that u inherits the regularity of /. I.e. there exists a constant K > such that, 
for any x' , z' £ R^ -1 and xn > 0, 

(6.14) - K\z'\ 2 < u{x' + z',x N ) + u(x' - z',x N ) - 2u(x' + z',x N ) < K\z'\ 2 . 
Then we repeat the 1-d proof essentially considering 

sup (u(x',x N ) - u(x',y N ) - C\x N - un^) ■ 

[0,+oo)x [0,+oo) 

Of course, a doubling of variables in x' is necessary to take care of the singularity 
of the measure, but using the W 2 '°° property in x' , we can go back to the 1-d 
computations without any difficulty. Let us just mention the key decomposition we 
use here. We rewrite the integrals with respect to first replacing the integrands 

by 

u{x' + z' , xn + zn) + u{x' — z' ,xn — zn) — 2u(x' , xn), 

and then by 

A 2 z ,u(x', xn + zn) + A 2 z ,u(x', xn — z N ) + 2A 2 N u{x' , xn), 

where 

2 z ,u(x' , xn) (u(x' + z' , xn) + u(x' + z' , xn) — 2u(x' , xn)^J , 

i(x' , xn) { u ( x ' ^ X N + z n) + u(x' , xn — zn) — 2u(x' , xn)^J ■ 

These expressions are not equal pointwise of course, but they give the same integrals 
because of the symmetry of /i*. We deal with the A 2 , -terms using (6.14), and the 
A 2 -term is treated as in the one dimensional case. Also note that we use a 
decomposition of VI into sets like I*" 1 x {zn ■ a < zn < b}, for a, b > 0, following 
the 1-d proof. 

Finally, concerning the nonsymetric part fj,#, we use as usual the fact that it is 
a controlled term since it is less singular. 



The existence is proved as in the proof of Theorem 6.2 □ 



Remark 6.1. The regularity results of the N = 1 and N > 2 cases are different. 
In the first case, the results is purely elliptic and we gain regularity. In the second 
case, the result is elliptic in the £jv-dircction while in the other directions we just 
use a preservation of regularity argument. It is an open problem to find an elliptic 
argument also in the ^'-directions. 

7. The limit as a -> 2" 

In this section we prove that all the Neumann models we consider converge to 
the same local Neumann problem as a — > 2 _ , provided that the nonlocal operators 
include the normalisation constant (2 — a). To be more precise, we consider the 
following problem 



(7.1) 



-(2 - a) J RN u a (x + rj(x, z)) - u a (x) dfi a + u a (x) = f(x) in fl , 

^ = in 90 . 

on 



where a £ (0, 2), r\ depends on the Neumann model we consider, and 

dfig _ g(z) 
dz ~ \z\ N + a 
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where g is nonnegative, continuous and bounded in WL N , g(0) > and g <G C 1 (B) 
for some ball B around 0. 

We prove below that the solution of (7.1) converge to the solution of the following 
local problem, 

!— aAu — b ■ Du + u = / in Q , 
— = mdil, 
on 

where 

IcJV-li I oiV-l I 

«:=ff(0)^^ and b := DgiO) 1 -^. 

In this section l^^ - 1 1 denotes the measure of the unit sphere in R N and Id^r the 
N x N identity matrix. 



Theorem 7.1. Assume (H^), i = 0...4 hold and let u a be the solutions of (7.1) 
for a € (0,2). Then, as a — > 2~ , u a converges locally uniformly to the unique 
solution u of (7.2). 

Before providing the proof, we introduce the following sequences of measures: 





= (2 




= (2 




-(2 


du tv 


= (2 



1 J \Z\N+" 



dz , 
dz . 



\A N 



(2 - a)n(y, z)iT)(y, z) dz , 



\N+a 



dz . 



where rj(y 1 z) i denotes the z-th component of the vector r)(y,z). Note that v\ and 
are matrix measures while v\ and v\ are vector measures. The localization 
phenomenon occuring as a — > 2 is reflected in the following lemma: 

Lemma 7.2. 

(a) As a — > 2~ , — a^old^r and v\ — bSo in the sense of measures. 

(b) For any sequence ak — > 2 andy^ — > x, there exist two vector functions d(x), b(x) € 
l" satisfying 

-a < a,i(x) < A and \b~i(x)\ < A for some A = A(<7, rf) < co , 
such that, at least along a subsequence, 

v a*m ~* diag(a(a;))(5 , v^ yk b(x)S , 
where diag(a(a;)) is the diagonal matrix with diagonal coefficients <ii(x). 
Proof. If S G (0, 1) is fixed, we notice first that, for any K > 1, 



0<(2-o) / M 2 f^£ < \\9\U5 

JS<\z\<K -! 



2-a _ K 2- a) ^ o as a 



so that the only possible limit in the sense of measure is supported in {0}. Similar 
calculations show that the same is true for all the measures v % , i — 2 . . . 4. 
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Coming back to v , we compute the inner integral as follows, 



(2 -a) / 

J\z\<6 

= 5 (0) (2 - a) [ z iZjT -^ + (2 - a) f 

J\z\<S \ z \ J\i 



'\z\<5 \ Z V 

The second integral vanishes as a — > 2 since 
f _g(z)-g(0) 

J\z\ 



\z\<b \ z \ 



(2 -a) 



>\z\<6 



% 3 



\N+a 



dz 



<C g {2-a) 



\N+a 



dz < C g (2-a 



6 3-a 

3 — a 



as a 



l\z\<6 \ z \ 

for C g = \\Dg\\ Lx ,( B/i y By symmetry, the first integral is zero for i ^ j, while for 
i=3, 



ff(0)(2-a) / 

J\z\ 



dz 



< s 1 \A N+a 



icJv-il r° r 



2+N-l 



dr 



I 



r Af+a 

a as a — > 2~ . 



This means that the measures {v^} concentrate to a delta mass i5o multiplied by 
the diagonal matrix aid at. 

Let us now consider the inner integral for each component of the measures v\ : 
using similar arguments, we have 

. g(z)-g{0) 



(2 -a) / 

JU 



\z\<6 
\z\<6 



\ z \N+a 



dz 



(2 -a)/ Zi (^9(0)) + (z,D 9 (z)-Dg(0)) dz 

J\z\ 



\N+a 



7 = 1 3 J \z 

(0)(2-a) / 

J\z 



dz 



dg 

dxi 



\<s \ z 
dz 

,\N+a 



1 l-IJV+a 



+ 05(1) 



+ 4 (1) 



- ^(0)4^ +0,(1) as a -2". 

Hence, by the definition of 6, ^ concentrates to b5o- 

We now come to the measures v z which is more complex to analyse due to the 
presence of the perturbation r](yk,z). We first notice that by using (H^), it follows 
that for i ^ j, 



\z\<8 



v(yk,z)iv(yk,z) 



g{z)dz 
3 \z\ N +° 



0. 



Then by (H*) \ij(y k ,z)\ < c v \z\, and we have 



0<g(0)(2-a) f r,{y k ,z) 

J\z\<8 



2 dz A 
iv+a 



< 5 (0) C 2(2-a) / \ z \z^- dz < g{ oy n \S 

J\z\<5 \ z \ 



2 I oJV-ll 
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So, the total mass of v 3 is bounded and, by the same arguments as above, it is clear 
that the support of v z shrinks to {0} (or the empty set). 
Then, we split the integral over \\z\ < 6} as follows 

(2-*)J ^,,«j^_/ w<i (...) + /„„ (■■■)-(*) + («). 

1 1 *N>-Vk,N z N<-Vk,N 

The first integral is easy to handle since T]i(yk, z) = z when zn > ~yk,N, 
(A \ (o n \ [ „2 9{z)dz 

z JV>-»fc,W 

2 g(z)dz 1 



z N >0 

The other integral has a sign and can take different values according to the structure 
of the jumps, but in all cases we see that the weak limit of v 3 can be written as 
a(x)5o where a(x) satisfies a/2 < di(x) < A. 

The measure v 4 is treated similarly: the total mass can be bounded by 

(2 -a) \r)(y,z)\ ,. n dz 

J\z\<5 \ z \ 

<c v C g (2-a) f \z\ 2 J J^ a =c n C g \S N - 1 \6 2 -", 

J\z\<5 \ z \ 

so that, up to a subsequence, there exists indeed a vector function b such that 
v a n y n ~> in the sense of measures, with ||6||oo < c r]Cg\S N ~ 1 \. The result then 
holds with A := \S N - 1 \c n max{C gi g(0)}. □ 

Remark 7.1. Note that in the censored case, a(x) = a/2 since the jumps below level 
— j/jv are censored, while a{x) = a by symmetry when the jumps are mirror reflected. 
Under our general hypotheses, different structures of the jumps (i.e. different 77's) 
lead to different a's which could in principle depend on x and the sequences ctk,yk- 
We will overcome this difficulty by using the extremal Pucci operator associated to 
a(x): for any symmetric N x N matrix A with eigenvalues (Aj) we define 

(7.3) M + (A):=^J2^ + A J2 

Ai<0 Ai>0 

Proposition 7.3. Let us define the half relaxed limits as a — > 2~ , 

u(x) := limsup u a (y) and u{x) := liminf u a {y) . 

a-t2,y^x a->2,y-tx 



Under the assumptions of Theorem 7.1 . u is a viscosity subsolution of (7.2), and u 
is a viscosity supersolution of (7.2). 

Proof. The proofs for u and u are similar, therefore we only provide it for u. We 
have to check that u satisfy the viscosity subsolution inequalities for the Neumann 
problem (7.2) at any point i6fl. There are two separate cases to check, (i) when 
x G CI and (ii) when x £ dfl. 

Step 1. Case (i) where x S fl, that is xm > 0. Let <j> be a smooth function and 
assume that x is a strict local maximum point of u ~ <j). By standard arguments 
there exists a sequence (y a ) a of local maximum points of u a — cf> such that y a — > x 
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as a — > 2 . Moreover, since xn > 0, by taking a close to 2, we can assume that 
Ua,N > $ for some small i5 > 0. By the subsolution inequality for u a at y a , 



-(2 - a) I 4>(y a + z)- 4>{y a ) - D<j>(y a ) ■ z d\i a - (2 - a) / D4>(lJa) ■ z d[i a 

J\z\<6 J\z\<b 

-(2 -a) / u a (P(y a ,z))-u a (y a )dfi a +u a (y a ) < f(y a ). 

J\z\>5 

We recall that the second integral of the left-hand side is well-defined : see the 
remark after Lemma 12. II 

We denote the three integral terms by I\, I2, and I3. Then 



h = -(2 -a) ( {{D 2 (j){y a ) + o s {l))z, z) dfi a 
J\z\<5 

= -(2 -a) (D 2 (/)(y a )z,z)dfi a + o s (l). 

J\z\<5 



Note that the o,s(l)-term is independent of a because the measure (2 — a)|z| 2 /z a 
has bounded mass. The symmetry of (x a implies that L <s ZiZj d[i a — and then, 
by Lemma [7. 2 1 we get 



h = -(2 - a)Tr(ZH0(y a )) / \z\ 2 d^ a + o s (l) = aA0(x) + o a (l) + o a (l) . 

J\z\<S 

Similarly we have 

h = -(2 -a) D<j)(y a ) ■ zd(i a = -(2- a)D(j)(y a ) zd[i a , 

J\z\<8 J\z\<S 

and by symmetry of (jL a and Lemma |7.2| we see that 



/„>= -(2- ( Ol.)o(//„)/ ^ ( M JV+ g i° ) dz = b- D<t>{x) + o a (l) + . 

'|z|<<5 l z l 



The 05(l)-terms are independent of a since the measures v\ of Lemma 7.2 have 
unformly bounded mass. For the last integral 7 3 , we use the boundedness of (u a ) a 
with respect to a to see that 

(7.4) \ h \<C ( 2- o) j M> ^<C^^ , s a^2. 

So we keep 6 > fixed and pass to the limit as a — > 2~~ (and y a — > x) to get 

-aA(j)(x) - b ■ D<j>(x) + u(x) < f(x) + o s (l) . 

Then, since S < xm could be arbitrarily small, we pass to the limit as 8 — > and 
get the viscosity subsolution condition for u at x. 

Step 2. Case (ii) where x € dfl, that is xn — 0. We again consider a smooth 
function cj> such that u — 4> has a strict local maximum point at x and, as above, we 
have a sequence (y a ) a of maximum points of u a — 4> such that j/ -}xasa4 2 _ . 
In this step we are going to prove that 

(7.5) mm M + (D 2 u(x)) - A\Du{x)\+u{x)- f(x) ; < , 

where M. + is defined in (7.3). We may assume that f^(a^) = — ~§^{ x ) > since 
otherwise (7.5| is already satisfied. Then for a close to 2, —-^^(y a ) > by the 
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2.1 



and this would contradict 



(2- 


a) 


(2- 


a) 


(2- 


a) 



continuity of D(j). We can also assume y a S f2, sinc e oth erwise y a S <90 and then 
§£(j/ a ) = -^{y a ) < for a close to 2 by Definition^ - 
our assumption. 

Therefore < y a ,N — > as a — > 2, and the subsolution inequality for u a takes 
the form 

4>(Va + V(y<x, z )) - HVa) - D(j)(y a ) ■ i](y a , z) d/i a 

\z\<5 

Dcj)(y a ) ■ r)(y a ,z)dii a 

\z\<5 

u a (P{y a ) - u a (y a ) d[i a + u a (y a ) < f{y a ). 

\z\>8 

We denote as before the three integral terms by I\, /2,^3- The compensator term 
I2 can be written as 

f dz 
h = -g{$)D<j){y a ) ■ (2 - a) / rj(y a , z) ^, N+a 



\z\<8 



z 



-DJ>(y a ).(2-a) [ ^z) 9 ^ ^ dz 

J\z\<5 \ z \ 

= l2.\ + ^2,2 ■ 

For symmetry reasons of both 77 and the measure, ^2,1 reduces to the scalar product 
of the 7V-th components, and it has a sign, 

d<j) . f dz 

i2,i = -ff(0)^ — {y a ){2-a) / r]{y a ,z) N > 0, 

OXN J\z\<8 \z\ N + a 



since <?(0), —g^{ya), and the 77-integral are nonnegative (see Lemma 2.2 (iii)). 
Thus we may drop the ^2,1 term from the inequality above and get that 

h + h,2 + h + u a {y a ) < f(y a ) . 

We now pass to the limit in this inequality as a — > 2 and hence y a ,N 0. The 
difference with Step 1 above, is that now y a converge to the boundary so that 



we cannot take a fixed < S < y a ,N as a — > 2. For the first integral, Lemma 7.2 



enables us to take subsequences «fc — >■ 2 and y a — > such that (dropping the 
subscript k for simplicity) 

g(z)dz 



h = -(2-a) / 4>{y a + rj(y a ,z)) - (j>{y a ) - D(f)(y a ) ■ rj(y a ,z) 

J\z\<8 
i,j J \z\<S 

= -J2 d lM x ) I d(isl v J l:3 (z)+o a (l)+os(l) 
= a i (x)dl l (f>(x) + o a (l) + 0,5(1) 

i 

> -M + {D 2 <t>(x)) + o a (l) + o s (l) . 

The last term I 3 can be treated as in Step 1 and vanishes as a — > 2. We are 
left with the I2 2 term and use again Lemma 7.2 this time for the measure ir. The 
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result is the existence of a vector b(x) such that along subsequences we have 
h,2 = D(j)(y a ) • (2 - a) 



\z\<8 \ z \ 

= D<j){x) ■ b{x) + o Q (l) > -A\D(f>(x)\ + o a (l) . 
Hence, passing to the limit a —¥ 2 in the above inequality, leads to 
-M + (D 2 (/)(x)) - A\D<f>(x)\ + u{x) - f(x) < 0, 



and (7.51 still holds. 



Step 3. We shall prove now that boundary condition ( 7.5 1 reduces to the condition 
< 0. Let us assume on the contrary that g^(ai) > for some point x at the 
boundary {xn = 0} and some smooth function <f) such that u — <f> has a maximum 
point at x. For any r, e > 0, we take a smooth, bounded function ip : R+ — > K+ 
such that 

ip(t) =r(t- 1 ^) for 0<t<e 2 /2. 
Since i/j(0) = and < ip for < t < e 2 /2, it follows that u(x) — 4>{x) — i/j(xn) has 



again a local maximum point at x. Hence (7.5) holds with 4>{x) + iP(xn) replacing 
(j>{x), i.e. 

-At; ®r(x)-r) <0, 
on / 



(7.6) 
where 



, a 2r 

mm(S(0) + -- 



E(<j>) := -M + (D 2 <j)(x)) - A\D4>(x)\ + u{x) - f(x) . 

Since we assumed that ^(x) > 0, we first fix r > small enough so that the 
inequality ^(x) — t > still holds. Then we can choose e > small enough to 



ensure that also 



E{4>) 



a 2r 
2^2 



At > 0. 



But then we contradict (7.6), and hence the boundary condition for u reduces to 



dip/dn < everywhere on the boundary. This concludes the proof of Proposi- 
tion El □ 



Proof of Theorem\7.1\ We have seen that u is a subsolution of (7.2) while u is a 



supersolution of the same problem. Since u < u on ft by definition and u > u on ft 
by the comparison principle for (7.2), we see that u = u on Q. Setting u := u = u 



on Q, it immediately follows that u is a continuous (since u is lsc and u is use) and 
the unique viscosity solution of (7.2). By classical arguments in the half-relaxed 
limit method, the sequence {u a ) a also converge locally uniformly to u. □ 



Appendix A. Blowup supersolution in censored case I. 

In this section we assume (H^j) and (H^) as in Section [5J Remeber that £1 := 
{(xi, . . . 7 xn) — (x' ,xn) '■ xn > 0}. First we show that in the censored fractional 
Laplace case (i.e. the censored alpha stable case), we can essentially take 

U{x) = — In xn 

as our blowup supersolution in assumption (U) in Section pjl 
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Lemma A.l. If dfi(z) = | z |w+c. for a e (0,1) and U{x) = — ln(xjv), then 

r dz 

-I[U\{x) = - U{x + z)-U{x) VJWTa >0 for x e Q. 

J x N +z N >a \ z \ 

Proof. We first change variables, z = -r-, to find that 

* JV 

Jx N +z N >0 v x/ \ z \ X N Jz N >-l \ z \ 

Now we are done if we can prove that 



r dz 
Jz N >-i M 
When JV = 1, we take 1 + z = e v and note that simple computations lead to 

e y dy f°° ^.m-a), u._. v 

|2sinh|| 1 +«' 



/ 

J — < 



V , y^La = I F{y)^ (1 - a) dy where F(y) 



Since F(y) is odd and 1 — a > 0, 

< -Ft-^e-^ 1 -") < F^e^ 1 "") for y > 0, 

and hence by symmetry J > 0. 

In the case N > 1 we introduce polar coordinates z — ry where r > and 
\y\ = 1, and we let dS(y) be the surface measure of the sphere \y\ — 1 in R N . We 
then find that 

~N-1, 



drdS(y) 
ln(l + ry N ) "Tf— ■ 



m/ / •/ r 

y J\v\=i,vn>o Jo J\v\=i,y N <o Jo 
The change of variables s = yjyr then leads to 

J=([ [°°+[ I 1 )sgn(y N )\y N rHl + S ) T ^dS(y) 

K J\y\=l,y N >oJo J\y\=l,y N <0J0 ' \ s \ 

f r°° ds 

= \y N \ a dS(y) / Mi + 8)^. 

J\y\=l,VN>0 J-l \ s \ 



The lemma now follows from the computations for we did for N = 1. □ 

We now generalize to a much larger class of integral operators with Levy measures 
fi such that d[i(z) ~ | z |w+ Q near \z\ = 0. In this case the blowup supersoution will 
be the modified log-function Ur defined as 

Ur{x) = Ur{xn) for x e ft, R> 1, 

where Ur is a (nonnegative) monotone decreasing C°°(0,oo) function such that 



U R {s) 



-ln(s) + |lni? if < s<R, 
if s > 2R. 



The main result in this appendix says that Ur will be the blowup "supersolution" 
of assumption (U) provided the Levy measure fi also satisfies: 
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([/)' For all R, e > there are r, c, K > and a £ (0, 1) such that 

(a) I ln(l + z N )[s a ^{sdz) - ° ^ > -e for s £ (0,r), 



l<2iV<-R 



(6) / ln(l + z N ) (i(sdz) > —K for s £ (r,R). 

J-l<z N <-± 

Theorem A. 2. Assume (H^) ; (H^)', and (U)' hold. Then the function Ur defined 
above satisfy the assumptions in (U). In particular, there is Rq > such that for 
any R > Rq there is Kr > such that 

-I[U R ){x)> -Kr in {x : < x N < R}. 

Before we prove this result, we show how assumption (U)' can be checked when 
fi is Levy measure whos restriction to {z : \zn\ < r} has a density 

(a 6(0,1), 

(A.i) f z =j$h where o^ei^n^;^), 

[\im z ^ g(z) = g(0) > 0. 

Note that the L 1 assumption makes g integrable near infinity and that L°°(M. N ) C 
L\M N ; 1+] f N+a ) for a > 0. 

Corollary A. 3. ///i /ias a density satisfying (A.I I, i/ien £/ie function Ur defined 
above satisfy the assumptions in (U). 



Proof. By Theorem A. 2 we have to check that (U)' holds. Part (b) follows from 
Holder's inequality since ln(l + s) £ 1, 0). Now we check part (a). Note that 

cdz g{sz) — c 
\z\ N + a = ~|7p 

Now choose c = g(0) and write 

cdz 



s ulsdz) — — = — dz. 

I - \N+a \~\N+a 



ln(l + z N )( s a ii(sdz) Vi 

l<ZjV<.R v \ z \ 

r dz 

>- sup \g(r)-g(0)\ |ln(l + zjv)| 



-s<r<_Rs ' ' ' ' ' ' J -Kz N <R,' ' "1^1*^+° 

Part (a) now follows since the last integral is finite for any R > 0, while the sup-term 
goes to zero as s — > by continuity of g at z = 0. □ 



Remark A.I. Assumption (A.I I also includes measures like 

Mi 

M = 

i=i 

where /Lti have densities satifying ( |A.1| for different (^ and 0;^. To see this, simply 
take a = max^ a, and g(z) = J2iLi 9i( z )\ z \ a ~ ai an( i n °te that g £ L 1 (M Ar ; 1+ |^+ a )■ 
We can even relax this assumption to include measures with zero or arbitrary neg- 
ative on provided that max^ a* remains in (0, 1). Finally we mention that we need 
some assumption to insure that /x does not give to much mass to the negative part 
of the integral —I[Ur\. In (A.I I we do this by requiring continuity at of g, but a 
carefull reader can extend this assumption to allow some discontinuities at 0. 
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Remark A. 2. In assumption (U) it is only the restriction of /u, to the set 

{z : -r < z N < Rr} n {z : -1 < z N < --} 

that plays any role. Hence if fj, satisfies (U)', by taking f small enough, so will fi+fl 
for any measure fx satisfying 

/ djl < oo and supp fi H {z : — 1 < zn < r} = for some f > 0. 

J\z\>0 

E.g. the delta-measure p, — Xa=i is °k if x at > 0. 



Proof of Theorem A. 2 First note that there is an Rq > such that 

dz 



Jr ■= / ln ( 1 + z Jv)H^+^ > °- 

^-1<zjv<-Ro l z l 

Indeed, in the proof of Lemma [A. 1[ we showed that = J > 0. The result then 
follows by the Dominated Convergence Theorem since the integrand is positive for 
z/v > and integrable. 

For any R > Rq, we note immediatly that Ur is a nonnegative decreasing func- 
tion which trivially satisfies the second part of (U) with lor(s) — Un(s) ' ^ e w ^ 
now check that Ur has the appropriate supersolution properties and hence com- 
plete the proof that Ur satisfies (U) under (£/)' . By the definition of Ur, we can 
write 

-I[U R ]{x) = f In (l + ^) fi(dz) + I R 

J -x N <z N <R,x N V X N / 

ln(l + y N ) n(x N dy) + I R , 



'-i<y N <R 

where Ir — — J z >Rx Ur{x + z) — Ur(x) fi(dz) > since Ur is decreasing. By 
assumption ([/)' we then find a r > such that for xm G (0, r), 

-x a N I[U R ]{x) >Jr+[ ln(l + (afr^tf dy) - > \j R > 0. 

J-Ky N <R V \y\ l+aJ 2 

When xtv € (r,R), another application of (U)' along with (H^)' leads to 
-I[Ur]{x) 

>([ +[ +[ )ln(l + — Wdz) 

v J-x N <z N <-^- J-^<2jv<i?.n^|<i J-^< ZJV <-Rn| 2 |>i ' v j; A r/ 

>-A- max |d* ln (l + s)| f max |i n (i + s )|/" d/i ( z ). 

Fivl 7| 2 |<i se(-|,f) J\ z \>i 

Since this last expression is bounded for xn € (r, R), this completes the proof. □ 

Appendix B. Estimates for the censored case II. 

Lemma B.l. Let(i(dz) = uprW, a £ (1,2), and define 9(x) = \x N f. If ft G (0,1) 
and iffi, iften 

f>0 if 0>a-l, 

l[0](x) — P.V. / <9(cc + z) - 0(x) fi(dz) 1 = i/ = a-l, 

J x N +z N >0 „ .» a 

I < ij p < a — 1. 
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Proof. First let (3 e (0, 1) and N — 1, and define 6{x) — \x\ f) . Note that the change 
of variables z — xz followed by 1 + z = e s reveals that 

l[0](x) = P.V. [ \x + z\ f3 -\x\ R " 

Jx+z>0 



= \xf- a P.V.f \l + zf-l 

Jz>-1 



1+a 



r°° 9sinh^ 
= \xf- a P.V. / - " * e^+^dx. 
1 1 J_ x |2sinhf|i+« 

When fj = a — 1, the integrand is odd and hence the integral is zero. For j3 > a — 1 
(/3 < a — 1) the exponential factor makes the integral positive (negative). Hence 
when /3 + 1 - a = 0, >0or<0, then I [6] = 0, > 0, or < respectively. 

When N > 1, a similar result holds for 0(x) = |ccjv|^. The idea is to work in 
polar coordinates. We set x = ry for r > and |y| = 1 and let dS(y) denote the 
surface area element of the iV-sphere \y\ = 1. We also use the change of variables 
ry N = fx N . 

I[§\(x) 

f dz 

= \x N + z N f -\x N f 

Jx N +z N >0 \ Z \ 

r^" 1 * dS(y) 



\xn + ry N f - \xn\ 

\y\ = l J x N +ry N >0 



i r+f pmo")*^ 

J\y\=l,y N >oJo J\y\=i, VN< o Jo J rl+a 

i»hl,!/»>0J0 J»=1,»«(<0J0 / v 7 l r l 



'\v\=i,yN>0 Jo J\y\=i,yN<0J0 

[J — a I I i- ■ -ifl ^ ^ 



= M p ~ a / M a dS(y)/ |l + f|"-l 

-'|y|=i,s/jv>o J-i I'M ' 

Here the first integral is just a positive constant while the second integral is the 
same we found in the N = 1 case. The conclusion is therefore as in that case: 
When /3 + 1 - a = 0, >0or<0, then I [6] = 0, > 0, or < respectively. 

□ 

Next we consider the two integrals 

B(a)= f \1 + zf - 1 dn(z), 

J-a-l 

G = j \l + zf + zf -2 dn(z), 

where a > 1, (3 € (0, 1), and d/j,(z) = jjfa for a E [1, 2). 
Proposition B.2. ///? = a — 1 i/ien i/iere is a k > sucft iftai 

5(a) + G < -k < 

for any a > 1. 

By continuity of the integrals in /3 we have the following corollary: 
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Corollary B.3. There is n > and ft > a — 1 such that 

B(a) + G< -n<0 

for any a > 1. 



To prove Proposition B.2 note that z + 1 < for z € (—a — 1, —a) (a > 1) and 



that the change of variable 1 + z = — e x in B(a) leads to 

r lna 9=:inVi^ a i / ,lna 9qinli^ 

* W |2co S hf|i+« e J / ln(Q _ 1) |2coshf|i+« 

For the G integral we have the following result. 
Lemma B.4. 

f ln2 2sinh^ „ M ,- , f 00 2sinh^ > 



_ ln2 |2sinhf|i+« 7 ln2 |2sinhf|i+° 



and if ft = a — 1, 

f°° 2 sinh 

G = 2 / 7— — ; — : — -77- ; d>X. 

J hl2 |2sinhf|i+« 
Proof. First note that by symmetry 

G = 2 lim / \l + zf -1 dfi{z). 

Then, since 1 + z > for z E (—1,1), the change of variable 1 + z = e x leads to 
G = 2 1im/ , 2sinh 3 , ef( 1+ ^"> dx. 



b^0+ 7(_ 0o ,l n 2)\(ln(l-6),ln(l+6)) |2sinh || 1+Q! 

Note that ln(l ± b) = ±6 + 0(b 2 ) and decompose the above integral as follows, 
/ (••■)<** 

J(-oo,ln2)\(ln(l-fc),ln(l+fc)) 

= (/ +/ "/ )(-)«** 

v «/(-oo,ln2)\(-6,4-6) J(-oo,ln2)\(ln(l-6),-h) ./(-cx>,ln2)\(ln(l+6),6) y 

Now since sinh a; = x + 0(x 3 ), the last two integrals are bounded by G6 2 p^ = 
Cb 2 ~ a for b -C 1, and we have 

G = 2 lim / {■■■) dx = 2\P.V. f + f )(•••) dx. 

*>->-0+ J(-oo,ln2)\(-6,6) \ J(-ln2,ln2) J (-00,- In 2) / 

A change of variables in the last integral then gives the first statement of the Lemma. 
The last part of the lemma follows since the integrand is odd when ft = a — 1, and 
hence the integral over (—In 2, In 2) vanishes. □ 



We also need the next lemma. 
Lemma B.5. Ifft = a-1, then 5(2) < - 
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Proof. We will show that 

f ln2 2 sinh ^ f ln2 2sinh ^ (3: + ln2) f°° G 

B(2) = / , * dx < / ^=4- — dx < / (• • • ) ds 



jo |2coBhf| 1 + Q -y |2sinh3^L2|i+« 7 " 2' 

The last inequality is trivial, and since sinh is an increasing function, the first 
inequaliy follows if we can show that 

x x ~\~ In 2 

cosh — > sinh for all x £ (0, ln2). 

2 ~ 2 v ' ' 

But this easily follows since fix") = cosh | — sinh x+ ^ n2 satisfy 

f'(x) = - sinh — cosh — — < for all x, 

J w 2 2 2 2 ~ 

/(ln2) = ^ 1 >0. 

□ 



Proof of Proposition B.2 Divide the integral B(a) into three parts 



(/•O /4n2Alna />ln2Vlna 

/ +/ +/ )(•••) dx. 

Jln(a-1)A0 Jln(a-l)V0 Jin 2 7 

Now we conclude since the first integral is negative, the second one is less than — ^ 



by Lemma B.5 and the last one is less than — § by definition of G. □ 
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